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INTRODUCTION 


The  AEROCRANE  is  a hybrid  heavy  lift  vehicle  which  shows  considerable 
potential  for  accomplishing  a variety  of  low  speed  heavy  lift  missions.  An 
extensive  experimental  and  analytical  study  of  the  trim  and  dynamic  stability 
characteristics  of  the  AEROCRANE  is  reported  in  Reference  1.  The  experi- 
mental investigations  described  in  Reference  1 were  designed  to  provide 
verification  of  equations  of  motion  to  predict  the  dynamic  stability  and 
control  characteristics  of  AEROCRANE  vehicles.  Excellent  agreement  between 
an  analytical  model  of  the  flight  behavior  and  the  experimentally  measured 
characteristics  of  a free-flight  dynamic  model  of  the  AEROCRANE  was  obtained. 

This  report  employs  the  analytical  model  of  Reference  1 to  study  various 
aspects  associated  with  the  precision  hover  capabilities  of  this  vehicle. 

The  vehicle  configuration  of  Reference  1 employed  conventional  heli- 
copter type  controls,  i.e.,  cyclic  and  collective  pitch.  It  has  been 
proposed  that  vertical  surfaces  or  winglets  could  be  mounted  on  the  tips  of 
the  rotating  wings  of  the  AEROCRANE  to  provide  effective  hover  control. 

This  report  compares  these  two  methods  of  control  in  hovering  flight  and 
in  addition  examines  other  questions  associated  with  precision  hover  of 
the  AEROCRANE,  including  a comparison  with  large  helicopters. 

First,  the  basic  features  of  cyclic  and  winglet  control  are  examined 
and  compared.  Then  using  a simplified  dynamic  model,  the  impact  of  these 
two  control  methods  on  the  precision  hovering  task  is  discussed. 

Since  the  basic  mission  of  an  AEROCRANE  vehicle  is  to  lift  and  move 
heavy  loads,  various  features  of  the  dynamic  motion  of  a sling  load  and 
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its  influence  on  the  vehicle  is  considered. 


The  hovering  equations  of  motion  for  the  AEROCRANE  with  a sling  load 
are  presented  in  nondimens ional  form,  such  that,  the  important  parameters 
are  evident  and  the  influence  of  vehicle  size  can  be  readily  determined. 

The  basic  hovering  dynamic  stability  characteristics  of  the  AEROCRAKS 
are  discussed  with  particular  reference  to  the  influence  of  vehicle  center 
of  gravity  position,  sling  load  attachment  point,  and  sling  load  length. 

Then  a series  of  simulation  studies  are  described  which  employed  an 
analog  computer  operating  in  real  time  to  examine  the  problem  of  hovering 
over  a spot  in  gusts  and  point  to  point  translation.  The  use  of  cyclic 
pitch  and  winglets  for  control  is  examined  and  the  results  compared  to 
results  of  a Heavy  Lift  Helicopter  simulation. 

Finally  the  hovering  control  and  stability  characteristics  of  the 
AEROCRAIfE  are  compared  with  conventional  helicopters. 
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CONTROL  METHODS 


Two  control  systems  are  examined  in  this  report  for  maneuvering 
the  AEROCRANE.  One  is  a cyclic  pitch  system  similar  to  that  used  on 
a conventional  helicopter  and  in  fact  was  the  system  employed  to 
control  the  dynamic  model  of  Reference  1.  A second  control  system 
involves  the  use  of  winglets  or  vertical  aerodynamic  surfaces 
mounted  on  the  tips  of  the  rotor  blades.  The  winglet  control 
characteristics  are  developed  in  Appendix  A.  The  cyclic  control 
primarily  produces  moments  on  the  vehicle  as  a result  of  a control 
deflection.  The  winglet  control  primarily  produces  a translational 
force.  It  should  be  noted  that  one  of  these  control  systems  cannot 
be  substituted  for  the  other.  That  is,  the  vehicle  can  be  flown  by 
a cyclic  system  only.  However,  it  does  not  appear  likely  that  the 
vehicle  could  be  controlled  by  the  winglet  system  only.  While  it 
is  expected  that  the  translational  force  control  would  be  effective 
in  maneuvering  the  vehicle,  it  would  not  be  capable  of  countering 
gusts  effectively  since  gust  inputs  will  produce  pitching  and  rolling 
moments  as  well  as  translational  forces.  In  this  section  the  nature 
of  the  control  forces  and  moments  are  discussed  and  then  a simplified 
model  of  the  dynamics  is  examined  in  order  to  obtain  some  insight  into 
the  relative  advantages  and  disadvantages  of  these  two  control  systems. 
A.  Cyclic  Pitch 

The  cyclic  system  involves  Imposing  a once  per  revolution  pitch 
variation  on  the  rotor  to  produce  primarily  a moment.  The  parameters 
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of  the  rotor  system  and  vehicle  are  such  that  the  rotor  blades  can  be 
assumed  to  be  rigid  (the  lowest  natural  frequency  of  the  rotor  blades  in 
flapping  is  considerably  larger  than  the  rotor  RFM)  so  that  a cyclic 
innut  with  its  maximum  and  minimum  along  the  fore  and  aft  axis  of  the 
vehicle  will  produce  a pitching  moment.  In  addition  to  the  moment 
produced,  a small  inplane  force  is  produced.  This  inplane  force  produces 
as  well  a moment  about  the  other  axis  of  the  vehicle  since  the 

center  of  gravity  of  the  vehicle  is  located  below  the  rotor  plane.  That 
is,  a longitudinal  cyclic  (Au)  input  produces  a pitching  moment  and 
a side  force.  The  side  force  acting  about  the  center  of  gravity  produces 
a rolling  moment.  In  order  to  examine  the  size  of  this  coupling  effect 
consider  the  theoretical  expressions  for  the  pitching  moment  and  side 
force  produced  by  longitudinal  cyclic  (Appendix  3) 
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If  the  center  of  gravity  is  located  a distance  r^  below  the  rotor  plane 
(the  center  of  buoyancy  is  assumed  to  lie  in  the  rotor  plane) 
then  the  ratio  of  the  rolling  acceleration  to  pitching  acceleration  produced 
by  longitudinal  cyclic  is 
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For  typical  AEROCRANE  parameters 
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and  for  a typical  r^  of  0.1  the  rolling  moment  produced  by  longitudinal 
cyclic  is  about  2 percent  of  the  pitching  moment. 

The  ratio  of  the  lateral  translational  acceleration  to  the  angular 
acceleration  in  pitch  is  given  by 
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Note  that  this  dimensional  quantity  depends  upon  the  size  of  the  vehicle 
as  it  is  directly  proportional  to  rotor  radius.  Non  dimensionalizing, 
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Taking  k^  = 1 

0 = 1.U2 

X = 0.43 

-/ 

— = 0.32  = - 0.C2 

2 

TMs  ratio  will  be  compared  with  the  winglet  characteristics  below. 

It  can  be  seen  that  this  ratio  will  always  be  small  since  it  is  pro- 
portional to  the  ratio  of  the  induced  velocity  to  the  tip  speed  (X^). 
Note  that  this  is  a side  force  in  conjunction  with  a pitching  moment. 
Owing  to  symmetry  it  can  also  be  seen  that  a rolling  moment  will 
produce  a horizontal  force  and  the  relative  magnitude  of  the  coupling 
terms  due  to  the  horizontal  force  arising  from  a lateral  cyclic  input, 
will  be  identical  to  the  side  force  terms  due  to  longitudinal  cyclic. 

B.  Winglet s 

The  winglet  is  thought  of  as  a control  that  primarily  produces 
a translational  force.  In  addition,  winglet  deflection  will  produce 
a moment  owing  to  the  fact  that  the  center  of  gravity  lies  below  the 
rotor  plane.  Appendix  A derives  expressions  for  the  winglet  control 
effectiveness.  The  winglet  control  is  azimuthally  uncoupled  in  the 
sense  that  a longitudinal  force  produces  only  a pitching  moment  and  no 
rolling  moment. 

The  ratio  of  the  translational  acceleration  to  the  pitching  acceler- 
ation can  be  obtained  from  the  expressions  given  in  Appendix  A as. 
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Again  it  may  be  noticed  that  the  dimensional  quantity  depends  upon 
the  size  of  the  vehicle.  Essentially,  if  a series  of  Froude  scaled  vehicles 
is  considered, as  the  vehicle  size  is  increased  the  translation  acceleration 
producec.  by  an  angular  control  deflection  is  invariant  with  size,  however, 
the  angular  acceleration  decreases  with  size. 

Non  dimens ionalizing  the  above  expression 


Placing  the  typical  values  used  above 
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For  a typical  r^  = 0.1 

t 

— = 25 

q 

Thus  the  winglet  control  emphasizes  the  translational  acceleration, 
producing  only  a relatively  small  angular  acceleration  and  essentially 
gives  an  effect  on  the  vehicle  which  is  the  inverse  of  the  cyclic 
control. 

One  further  comparison  is  of  interest.  That  is,  the  question  of 
how  much  cyclic  is  required  to  cancel  the  angular  acceleration  of  the 
winglet  such  that  the  winglet  control  produces  translational  acceleration 
only.  The  angular  acceleration  produced  per  unit  cyclic  is 
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The  angular  acceleration  produced  by  the  winglets  is  given  in  Appendix  A 
as 

1 SMsw 


36 


A6  = (i  p(QR)^  Sm  r a^)  A6 
T-/  00  c ' / ^ w o ( 

I c I 


Taking  the  ratio  of  these  t-vo  terns 
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Introducing  the  definition  of  a = — r,  and  assuming  that  the  winglet 

TTK 

is  of  rectangular  planform  with  a chord  equal  to  the  blade  chord,  and 
taking  the  number  of  blades  as  four,  this  ratio  can  be  expressed  as 


— = I ""n  X (— ) 


A6 

(—  ) 
X A A > 


o ''  ' a ' Mi  'AAi3 


Inserting  numbers  into  this  expression,  assuming  that  b,^  = 0.2^(a  value 
shown  in  Appendix  A as  suitable  for  trimming  the  AEROCRANE  in  forward 
flight)  , and 


— = 0.8 

a 


(to  account  for  the  finite  aspect  ratio  of  the  winglet)  and  using  the 
2C  „ 

value  of  = 0.062  corresponding  to  a C t = O.OOSh  from  Appendix  B: 
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Thus  for  a typical  r^  = 0.1,  the  angular  acceleration  produced 
by  winglet  deflection  is  roughly  7 percent  of  the  angular  acceleration 
produced  by  cyclic  deflection. 

The  total  angular  acceleration  produced  by  the  winglet  can  be  ex- 
pressed as 


3-  = l+  -S.  = i + 


AAi3 

1 

0-7  A6  ^ 


It  is  possible  to  eliminate  the  angular  acceleration  produced  by 
the  winglet  control  by  mixing  cyclic  into  the  winglet  control  stick. 

The  amount  required  is  obtained  by  taking 

AAi3  = A6h 

Again  for  a typical  center  of  gravity  location  r^  = 0.1,  the  coupling 

required  is 

K = - 0.07  — — 

**  °winglet 

indicating  that  the  cyclic  control  produces  large  moments  compared  to 
the  winglet.  This  gearing  would  of  course  also  result  in  a translational 
force  along  the  other  axis,  that  is,  use  of  a coupling  of  longitudinal 
cyclic  to  eliminate  the  pitching  moment  produced  by  cosine  or  transla- 
tional force  winglet  deflection  will  result  in  a side  force. 

Above  we  have  developed  expressions  for  the  ratio  of  translational 
acceleration  to  angular  acceleration  of  the  two  controls.  For  typical 
parameters , for  cyclic  control 
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For  winglet  control 
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Now  with  the  cyclic  control  coupled  to  the  winglet  control  to 

eliminate  the  angular  acceleration  produced  by  the  winglet,  it  is  implied 

/ 

that  = - q^'  and  therefore  the  ratio  of  lateral  acceleration  produced 
by  cyclic  deflection  to  translational  acceleration  produced  by  winglet 
deflection  is 

t 

V 

— = - 1.28  X r 

— / so 

u„ 

Again  inserting  typical  values  of  r^  = 0.1  and  X^  = - 0.06 

/ 

V 

— = 0.0077 

u; 

Thus,  while  mixing  cyclic  with  winglet  deflection  to  eliminate  the 
angular  acceleration  produced  by  winglet  deflection  does  produce  a 
control  coupling  in  the  sense  that  a lateral  force  is  produced  in  addition 
to  the  longitudinal  force  arising  from  winglet  deflection,  the  effect 
is  essentially  negligible.  In  effect,  the  above  ratio  is  equivalent  to 
shifting  the  phase  of  the  winglet  control  by  0. 0077  radians  (.4U  degrees) 
to  obtain  a translational  force  only  with  the  cyclic  mixing  to  eliminate 
the  pitching  moment.  Precisely  speaking,  there  is  still  a small 


10 


moment  remaining  about  the  other  axis  (roll  in  the  case'''^f  this 
example)  owing  to  the  fact  that  the  cyclic  has  been  mixed  with  the 
winglet  deflection,  however,  this  effect  should  be  very  small. 

C.  Summary 

Thus  to  summarize  the  results  of  this  section,  a cyclic  control 
input  primarily  produces  a moment  on  the  vehicle  and  secondarily  pro- 
duces a small  translational  force.  A winglet  control  primarily  produces 
a translational  force  and  secondarily  produces  a small  moment.  For 
typical  aEROCRAHE  parameters  only  a small  level  of  cyclic  mixing  with 
winglet  deflection  is  required  (less  than  0.1  deg  cyclic  per  deg  winglet) 
to  eliminate  the  moment  produced  by  winglet  deflection.  This  ratio  is 
of  course  dependent  upon  the  center  of  gravity  location  on  the  vehicle 
and  would  be  of  the  order  of  1 if  the  center  of  gravity  of  the  vehicle 
was  located  on  the  gondola  as  a limiting  case. 


11 


HOVERING  CONTROL  CONSIDERATIONS 


In  order  to  obtain  some  insight  into  the  basic  mechanics  of  hovering 
control  of  the  AEROCRANE  a simplified  model  of  the  hovering  dynamics  is 
examined.  In  this  model,  to  bring  out  the  central  features  of  hovering 
control,  the  long  time  dynamics  and  the  sling  load  dynamics  are  neglected. 

That  is,  it  has  been  shown  in  Reference  1 that  the  basic  modes  of  motion 
of  the  AEROCRANE  include  a rapid  angular  mode  whose  characteristics 
depend  primarily  on  the  rotor  angular  damping  and  the  gyroscopic  effects 
of  vehicle  and  rotor  rotation.  In  addition,  two  long  period  oscillatory 
modes  exist  with  frequencies  primarily  determined  by  the  velocity  stability 
(M^),  dihedral  effect  (L^)  and  the  aerodynamic  coupling  terms  and  L^. 
Typical  periods  for  these  motions  as  shown  in  a later  section  of  this  report, 
are  of  the  order  of  ^ seconds  and  ^ seconds.  For  a large  AEROCRANE  with 

'^4  .2 

a rotor  angular  velocity  of  the  order  of  one  radian  per  second,  the  cor- 
responding periods  of  these  slow  motions  are  60  and  40  seconds.  Consequently, 
it  would  not  be  expected  that  these  longer  period  modes  would  be  particularly 
significant  in  precision  hover  control.  Thus,  in  the  following  analysis  the 
speed  dependent  derivatives  are  neglected  with  respect  to  their  influence  on 
the  vehicle  dynamics.  Only  the  translational  gust  response  can  be  examined 
with  this  simplified  model.  The  rapid  angular  response  mode,  prima. ily 
determined  by  the  terms  andcuj,  is  retained.  This  mode,  characterized 
by  a period  of  the  order  or  5 seconds  for  Q = 1.0  rad/sec.,  is  the  one 
primarily  involved  in  maneuvering  near  hover. 

Thus  the  following  equations  of  motion  are  employed  in  this  section 


to  examine  the  hovering  dynamics  and  control  of  the  AEROCRAKE. 


T 

where  T]  = — . The  symmetry  of  the  vehicle  has  been  employed  to  express  the 
m' 

lateral  derivatives  in  terms  of  the  longitudinal  derivatives.  Note  that  the 
derivative  is  entirely  due  to  the  winglets.  That  is,  the  centerbody  drag 
is  assumed  to  be  proportional  to  the  square  of  the  velocity  and  therefore  does 
not  contribute  to  the  linearized  model  when  hovering  is  the  equilibrium  con- 
dition. It  was  shown  in  a previous  section  that  the  translational  forces  due 
to  cyclic  are  small  and  therefore  the  terms  are  neglected. 

The  characteristic  equation  for  the  dynamic  motion  of  the  vehicle  described 
by  this  simplified  model  is 

= sMs  - [(s  - + u)/] 

The  short  term  dynamics  are  described  by  the  quadratic  factor  in  the  square 
brackets  and  the  long  term  dynamics  are  given  by  the  first  two  factors  which 
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correspond  to  two  roots  at  the  origin  and  two  real  roots  at  X . Without 

u 

winglets,  = 0 and  the  long  term  approximation  to  the  dynamics  consists 
of  four  roots  at  the  origin. 

The  transfer  functions  describing  the  translational  response  of  the 
vehicle  can  be  written  as 

(s)(s  - X^)  [s  AX  [ (s  - + u)5^}  - AMTl  (s  - M^)] 

A 

c 

(3)(s  - X^)  [-  AM  uja  Tl] 

^ 

c 

where  the  influence  of  various  inputs  is  shown  by  the  following  expressions: 

AX  = X.  6 + X Ua 

5 c u 
c 

AM  = M,  Ais  + M.  6 

.IS  Oc 

Only  longitudinal  inputs  have  been  included  at  this  point.  Owing  to  symmetry 
the  lateral  response  problem  will  be  similar  so  only  the  longitudinal 
is  discussed. 

First  we  contrast  the  basic  features  of  winglet  control  to  cyclic 
control  for  hovering  position  control.  In  the  first  instance  it  is  assumed 
that  for  the  winglet  control  case  the  cyclic  is  coupled  to  the  winglet 
such  that  no  pitching  moment  is  applied  to  the  vehicle  as  a result  of  a 
winglet  input.  That  is,  the  cyclic  is  geared  to  winglet  motion  such  that 

M»  Ai,  + Ms  6 =0 

resulting  in  a gearing 
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This  is  a lew  level  of  gearing  as  noted  in  an  earlier  section. 
The  longitudinal  velocity  response  to  winglet  input  becomes , 
after  inserting  the  expression  for  the  characteristic  equation  into 
the  transfer  function 


Consequently  the  expression  for  longitudinal  position  response  is 


If  it  is  assumed  that  the  pilot  attempts  to  control  longitudinal  position 
by  moving  the  winglet  control  proportional  to  longitudinal  displacement, 
a closed  loop  system  which  is  stable  with  a closed  loop  frequency  dependent 
upon  gain  is  obtained  as  shown  in  Figure  la. 

Before  considering  this  situation  in  more  detail,  consider  the 
contrast  with  cyclic  control  of  position.  From  the  transfer  functions 
given  above  the  longitudinal  velocity  response  to  cyclic  is  given  by 

- M,  T1  (s  - M ) 

u ^ ^ 

s=  [(s  - 

has  been  set  equal  to  zero  since  it  is  assumed  that  the  vehicle  is 
not  equipped  with  winglets.  The  transfer  function  for  longitudinal 
position  is 

- M,  Tl  (s  - M^) 

3^  ^ If ^ 

sM(s  - 
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Again  assuming  that  the  pilot  attempts  to  control  longitudinal  position 
with  cyclic,  the  closed  loop  dynamics  shown  in  Figure  lb  result. 

The  closed  loop  system  is  unstable,  illustrating  the  classical  problem 
of  controlling  the  translational  displacement  of  a hovering  vehicle. 

That  is,  the  pilot  must  supply  a considerable  amount  of  lead  in  order 
to  stablize  the  vehicle.  The  feedback  transfer  function  required  to 
to  stabilize  the  vehicle  must  be  of  the  form 

Ais  = - Kp  (t^i  s + 1)  (t,_2  s + 1)  X 
With  sufficient  gain  and  a double  lead  this  closed  loop  system  is 
stable  as  shown  in  Figure  Ic.  The  control  of  translational  displacement 
with  cyclic  pitch  is  a difficult  task. 

The  double  lead  term  can  be  interpreted  in  another  way  by  noting 
the  relationship  between  translational  displacement  and  vehicle  attitude 
from  the  equations  of  motion.  The  translational  force  equation  gives 
the  following  relationship  between  vehicle  attitude  and  displacement. 

Tie  = - s®  X 

Thus  the  double  lead  term  may  be  interpreted  as  a requirement  that  the 
pilot  provides  an  attitude  feedback  or  inner  loop  plus  an  outer  loop 
involving  translational  displacement  and  velocity  feedback. 

The  interest  in  the  use  of  winglets  relates  to  hovering  control 
where,  by  direct  control  of  translational  force,  two  integrations  are 
eliminated  from  the  transfer  function  relating  longitudinal  displacement 
to  control  action  simplifying  the  pilot's  control  task. 
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It  should  also  be  noted  that  if  the  winglet  configuration  does 
produce  an  appreciable  pitching  moment  in  addition  to  translational 
force  then  the  winglet  control  problem  becomes  somewhat  more  complex 
as  shown  below.  The  transfer  function  for  longitudinal  displacement 
including  the  influence  of  winglet  pitching  moment  is 

Xg  [s  {(s  - M (s  - M )] 

X _ c 

(s^)  (s  - X^)  [{s  - ] 

Figure  Id  shows  the  effect  of  the  ratio  - — on  the  location  of  the 

c 

zeros  of  this  transfer  function.  In  the  limit  of  a very  large  ’/alue 

of  the  ratio  - — the  zeros  approach  the  cyclic  case.  Now  consider 

the  effect  of  translational  position  feedback  with  some  significant 

value  of  the  ratio  - — as  shown  in  Figure  la.  Thus,  if  the  winglet 

«c 

» 

produces  an  appreciable  pitching  moment  as  well  as  a translational  force, 
the  hovering  control  of  position  becomes  more  complex  and  the  closed 
loop  dynamics  are  marginally  stable  unless  the  pilot  provides  a single 
lead  as  shown  in  Figure  If.  If  the  vehicle  configuration  is  such  that 
a winglet  input  produces  a significant  pitching  moment,  then  it  is  de- 
sirable to  provide  cyclic  mixing  such  that  a winglet  input  provides 
translational  force  only. 
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This  simple  model  is  now  employed  to  examine  the  translational 
response  of  the  vehicle  to  gusts.  Owing  to  the  fact  that  the  trans- 
lational velocity  pitching  moment  derivatives  have  been  neglected,  the 
model  is  valid  for  translational  investigations  only.  The  gust  response 
of  the  complete  dynamic  model  is  described  in  a later  section.  Thus, 
the  transfer  function  relating  translational  displacement  to  gust 
velocity  is. 


and  Xg  are  simply  related  as  shown  in  Appendix  A since  both 
c 

terms  arise  from  the  winglets. 

Let 


and  noting  that 

c 


Y -i  — =02t5-' 

\ ~ , b(,  ^ ^ e 

c m c 


and 

Xu  = -ne' 


The  translatici-u-  ^ust  response  transfer  function  can  be  written  as 


X 

Ua 


1 

^s  + e s + e K 
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where  the  Laplace  variable  has  been  nondimens ionalized  by  the  rotor 
angluar  velocity 


and  K is  a nondimens ional  gain, 

K = R K 

X X 

is  the  dimensional  feedback  gain  in  radians  per  foot. 

The  mean  square  value  of  the  displacement  response  for  a typical 
gust  input  spectrum  is  now  evaluated.  A reasonable  model  for  the 
power  spectrum  of  horizontal  gusts  at  low  wind  velocities  can  be 
obtained  by  passing  white  noise  through  a low  pass  filter.  The 
power  spectrum  of  a gust  with  mean  square  gust  velocity  u^^  is  given 


by 


ujcH 


s*  2 


uu  X2 


1 + C— ) 

(JU 

c 


where  ou^  is  the  cutoff  frequency  of  the  low  pass  filter  and  is  usually 
taken  to  be  0.314  rad/sec  as  a reasonable  match  to  experimental  data 
on  gust  spectra. 

The  power  spectrum  of  the  translational  displacement  can  then  be 
calculated  from  the  transfer  function  given  above. 

Gx(uj)  = 1;^ 

The  mean  square  value  of  the  displacement  of  the  vehicle ,x^  is  given  by 
^^  = 11  G^(tu)  dm 
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Explicitly 


SV 


K ou 
X c 


1 * (- 


e ' + uj 


or 


The  function,  f,  is  plotted  in  Figure  2,  for  a typical  value  of  s' . It 
can  be  seen  that  the  critical  parameter  is  the  dimensionless  feedback 
gain,  K^,  and  the  mean  square  displacement  is  not  particularly  sensitive 
to  (ju^.  Thus,  the  influence  of  size  of  the  vehicle  enters  primarily 
through  the  rotor  RIM.  For  the  dynamic  model  of  Reference  1,  with 
Q = 3-1^  rad/sec  and  a root  mean  square  gust  level  of  5 fps , the  value  of 
f = 0.39  yields  a root  mean  square  (RMS)  displacement  of  1 ft.  Figure  2 
indicates  that  a dimensionless  gain  of  the  order  of  2 is  required  for 
this  RMS  displacement.  To  determine  the  dimensional  gain, 

K 

and  for  the  model  of  Reference  1,  R = 19.64  ft.  and  therefore  the  dimen- 
sional gain  is 

K°  = 5.8  degrees  winglet  deflection  per  foot  displacement. 
This  gain  corresponds  to  a closed  loop  natural  frequency 


R e'  K = 2.i+2  rad/sec. 

X 


or  a 2.6  second  period.  If  the  l6  ton  vehicle  is  examined,  with 
0 = l.lU  rad/sec,  for  a root  mean  square  displacement  of  1 foot, 
f = 0.046  is  required  and  this  corresponds  to  a dimensionless  gain 
= 6.  The  dimensional  gain  (R  = 136.5  ft.)  required  is  about  one- 
half  the  value  for  the  smaller  vehicle.  That  is  K°  = 2.51  degrees  per 
foot.  This  example  illustrates,  as  indicated  also  in  Reference  4,  that 
in  order  to  maintain  position  in  gusts  a relatively  high  gain  displace- 
ment feedback  is  required. 

For  the  vehicle  without  winglets,  the  linearized  model  would 
indicate  no  response  to  gusts  since  the  centerbody  forces  depend  upon 
the  square  of  the  velocity.  The  analog  computer  model  included  this 
nonlinear  effect  and  the  complete  gust  response  of  the  vehicle  based  on 
the  computer  model  is  discussed  in  a later  section. 
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a.)  Winglet  Control  of  Displacement  (No  Winglet  Pitching  Moment) 

Control  Law:  6 = K x 

c w 


b.)  Cyclic  Control  of  Displacement 

Control  Law:  Ai.  = - K x 

c 


Figure  1:  Comparison  of  Winglet  and  Cyclic  Control  in 
Hovering  Flight.  Simplified  Dynamic  Model. 
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e.)  Winglet  Control  of  Displacement  Including  Moment  Effect  of  Winglets 
Control  Law:  6 = - K x 


f.)  Winglet  Control  of  Displacement  Including  Moment  Effect  of  Winglets 
and  Pilot  Lead 

Control  Law:  6 = - K (x,  s + 1)  x 

c c ^ 


Figure  1:  Continued. 


SLING  LOAD  DYTJAMICS 


In  this  section  seme  basic  features  of  the  sling  load  motion  are 
examined.  Of  particular  interest  is  the  influence  of  the  presence  of 
a buoyant  force  in  addition  to  aerodynamic  thrust.  Consequently  ’/arious 
results  obtained  are  contrasted  with  a hovering  helicopter  for  which 
the  buoyant  force  equals  zero  (3  = O).  The  analj'-sis  presented  is  based 
on  a simplified  model,  such  that,  it  not  only  applies  to  the  AEROCRANE 
but  to  any  hovering  vehicle  which  obtains  part  of  its  lifting  force 
from  buoyancy.  Of  particular  interest  therefore  is  the  influence  of 
the  buoyant  force  ratio,  0,  and  the  sling  load  mass  ratio,  u , on  the 
sling  load  and  vehicle  dynamics. 

It  is  shown  in  a later  section  of  the  report  that  the  sling  load 
geometry  (the  location  of  the  sling  load  attachment  point,  and  the 
sling  length)  has  only  a small  influence  on  the  dynamic  stability  character- 
istics of  the  vehicle  (i.e.,  the  characteristic  roots),  that  is,  in  a 
dynamic  sense,  the  sling  load  is  not  strongly  coupled  to  the  vehicle. 
Therefore,  it  is  possible  to  examine  the  dominant  effects  of  the  sling 
load  using  an  approximate  model. 

It  is  assumed  that  the  sling  load  motion  and  its  effect  on  vehicle 

motion  can  be  estimated  by  considering  the  system  in  the  following  way. 

The  interaction  between  sling  load  motion  and  vehicle  motion  is  approxi- 
mated as  a forced  response  problem.  That  is,  the  pitch  attitude  of  the 
vehicle  is  assumed  to  be  the  input  and  the  sling  load  motion  and  the 

vehicle  translation  motion  are  assumed  to  be  outputs.  Thus,  the  motion  of 

the  complete  dynamic  system  is  reduced  to  an  examination  of  a relatively 
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simple  frequency  response  problem.  Since  vehicle  pitch  attitude  is  taken 
as  the  input,  it  has  been  assumed  that  cyclic  pitch  is  being  employed  for 
control  and  the  input  frequency  characterizes  the  pilot's  control  action. 

It  has  been  shown  by  frequency  response  experiments  with  the  complete 
dynamic  model  of  the  AERCCRANE  that  this  simplified  approach  is  valid. 

With  this  simplified  model,  the  following  discussion  applies  to  any 
heavy  lift  vehicle  which  is  controlled  by  moments  and  derives  an  appreciable 
portion  of  its  lift  from  buoyancy,  since  only  the  translational  force  and 
sling  load  equations  are  involved  in  the  analysis.  If  is  neglected, 
the  horizontal  force  equation  as  given  in  Appendix  C is 

m'  = - (w^  - Fs  ) 9 + Ws , e, 

The  sling  load  equation  is 

- Ug  - g (0  + 9^  ) - 0 (Z,  + Z J - ©L  Z^  = 0 

The  horizontal  acceleration  is  eliminated  from  the  sling  load  equation,  using 
the  translational  force  equation  to  obtain  an  equation  relating  load  deflection 
relative  to  the  vehicle  (9^)  and  vehicle  attitude  (9).  The  geometry  is  shov/n 
in  Figiure  3-  The  resulting  equation  is. 


Nondimensionalizing  the  time  by  the  sling  load  frequency  and  introducing  the 
dimensionless  parameters  u and  0 describing  the  sling  load  weight  ratio  and 
the  buoyancy  ratio,  the  above  equation  becomes 


This  equation  can  now  be  viewed  as  describing  the  response  of  the  sling 
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load  to  vehicle  attitude  changes.  Taking  the  Laplace  transform  and  reorganizing 

Za 

I (1  ^ ^ 

9, 


T 


^ a 

(1  + s"  + — ^ 

(1  + I) 


I (1  + 1 + : 


(1  + 1) 


Consider  first  for  reference  a conventional  helicopter  with  0=0 


s®  + ( 1+u, ) 

If  the  load  is  attached  at  the  center  of  gravity,  such  that  = 0 

( 

|s^+(l+u) 

Now  consider  the  frequency  response  of  9^  to  9 inputs,  i.e.,  let 
s = icu  recalling  that  the  time  has  been  nondimensionalized  by  the 
sling  load  frequency.  It  is  clearer  at  this  point  to  examine  the 
sling  load  motion  with  respect  to  space.  Let 

9la  = 9l  +9 


The  transfer  function  for  the  absolute  load  motion  is  therefore 


9 


1+u 

s®+  (l  + M.) 
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With  s = iu) 


21  = (1  + u) 

® -35^+(l+lj,) 

at  low  frequencies  as  u3  0, 

9la 

indicating  that  the  motion  of  the  load  in  space  is  the  same  as  the 
helicopter  motion  and  at  high  frequencies  the  absolute  load  motion 
tends  towards  zero,  that  is,  the  load  remains  fixed  in  space.  This 
frequency  response  is  shown  in  Figure  U. 

Now  consider  the  case  in  which  the  load  is  attached  some  distance 
below  the  center  of  gravity.  The  transfer  function  is 

Qla  ■ ^ " + (1  +M.) 

s ® + (l  + U ) 


The  low  frequency  behavior  is  unchanged,  however  there  is  a difference 
in  the  high  frequency  response 


0L* 

9 

UJ  -• 


00 


Za 


and  the  complete  frequency  response  is  given  by 


Za 

— i ^ + (1  + U ) 

^ L 

= + (1  + u) 
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Thus  for  the  nose  up  motion  of  the  helicopter,  the  load  moves 

aft  of  the  vertical  by  an  amount  dependent  upon  the  distance  between 

the  center  of  gravity  and  the  load  attachment  ooint.  Note  that  in  the 
Za 

case  where  — equals  one,  although  the  amolitude  of  the  load  is  essentially 
Zl 

independent  of  frequency,  the  phase  shifts  by  l80°  so  that  the  load  motion 
is  in  fact  different  for  forcing  above  and  below  the  natural  frequency  of 
the  load. 

Now  consider  the  case  with  buoyancy.  The  transfer  function  for  the 


absolute  load  motion  in  response  to  vehicle  attitude  is 


(1 


The  high  frequency  response  is  unchanged,  however,  the  low  frequency 
response  becomes 


U)  -•  o 
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It  can  ba  seen  that  in  the  fullybuoyant  case  (3  = 1 + ij, ) the  load 
remains  vertical,  that  is,  the  absolute  load  motion  tends  towards  zero. 

Thus,  the  effect  of  adding  a buoyant  force  to  a hovering  vehicle  is  to 
tend  to  reduce  the  low  frequency  response  of  the  load  motion.  This  is, 
of  course,  a consequence  of  the  fact  that  with  a buoyant  force  there  is 
less  translational  acceleration  as  a result  of  tilting.  Figure  5,  shows 
the  frequency  response  of  the  load  in  this  case . 

To  summarize,  the  high  frequency  (above  the  sling  load  natural 
frequency)  maneuvering  excites  the  motion  of  the  load.  If  the  load  is 
attached  at  the  center  of  gravity  of  the  vehicle  the  load  tends  to  remain 
vertical  and  with  the  load  below  the  center  of  gravity  the  absolute  motion 
of  the  load,  with  respect  to  space,  is  increased.  Thus,  rapid  maneuvering 
will  tend  to  excite  the  load  motion  which  has  little  damping  and  will  tend 
to  continue  to  oscillate.  In  the  slower  maneuvering  of  the  vehicle  with 
a buoyant  force  present,  the  load  tends  to  remain  near  the  vertical.  It 
should  also  be  noted  that  the  use  of  winglets  will  tend  to  excite  the 
load  as  winglet  deflection  will  produce  a horizontal  force  and  consequently 
the  load  will  experience  an  acceleration  and  commence  swinging. 

Another  important  aspect  of  the  influence  of  the  sling  load  on  the 
response  of  the  aEROCRANE  can  be  seen  by  examining  the  translational 
acceleration  response  to  vehicle  attitude.  This  quantity  can  be  thought 
of  as  a measure  of  the  translational  response  of  the  vehicle  to  cyclic  pitch 
with  a tight  attitude  loop  such  that  it  is  reasonable  to  assume  that  cyclic 
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pitch  produces  attitude.  The  horizontal  acceleration  response  to  attitude 
can  then  be  examined  using  only  the  sling  load  equation  and  the  translational 
force  equation,  that  is,  pitch  attitude  is  considered  to  be  an  input.  It 
should  be  noted  that  most  of  the  features  of  this  response  discussed  below 
are  characteristic  of  the  response  of  any  heavy  lift  vehicle  which  obtains 
an  appreciable  amount  of  its  lift  from  buoyancy. 

The  translational  force  equation  as  given  above  is 


u 


s 


The  important  point  to  note  about  this  equation  is  the  fact  that  when  there  is 
an  appreciable  buoyant  force  Fg,  the  influence  of  vehicle  attitude  (9)  upon  the 
horizontal  acceleration  becomes  small  compared  to  the  influence  of  sling 
load  deflection  (9^  ).  Expressing  this  eq’uation  in  terms  of  dimensionless 
ratios , 


s 

g 


-I 


For  the  AEROCRAUE  parameters  assumed  here, 

0 = 1.42 

u = 0.84 
u 

— = + 0.24  9 + 0.49  0^ 
g 

The  effect  of  sling  load  deflection  is  about  twice  the  influence  of  vehicle 
attitude  and  in  addition,  the  effect  of  vehicle  attitude  has  changed  sign 
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compared  to  a helicopter.  Thus,  the  translational  acceleration  response 
to  vehicle  attitude  will  depend  to  a large  extent  on  the  sling  load 
response  to  attitude  discussed  above.  Consider  this  equation  for  a 
helicopter  where  0=0.  is  taken  to  be  0.75  for  a typical  heavy- 
lift  helicopter,  the  corresponding  equation  is, 

’ig 

— = - 9 + 0.75  9, 
g 


The  attitude  term  changes  sign  and  is  larger  than  the  sling  load  de- 
flection term.  In  order  to  understand  further  the  implications  of 
the  above  results  consider  a low  frequency  motion  of  the  vehicle,  in 
which  the  angular  acceleration  terms  can  be  neglected  in  the  sling 
load  equation.  In  this  case,  the  sling  load  equation  can  be  solved 

for  sling  load  amplitude, 

u 


9l 


Substituting,  for  the  AEROCRAUE 
u 

— = - 0.167  9 

g 


For  the  helicopter  with  no  buoyant  lift. 


g 


Thus,  for  slow  maneuvering  of  the  buoyant  vehicle  with  a sling  load 
only  about  ^ the  acceleration  is  obtained  as  compared  to  the  helicopter. 
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Now  the  more  general  case  is  considered,  that  is,  an  expression  is 
derived  for  the  frequency  response  of  translational  acceleration  to  vehicle 
attitude.  The  relationship  between  the  sling  load  motion  with  respect  to 
space  and  the  vehicle  attitude  was  developed  above  as 


e 


_ 1:  ~3  ,(1  - P +a) 

(1  + I) 


+ 


(1  I 

(1  +|) 


This  expression  can  be  used  with  the  sling  load  equation  to  eliminate  the 
sling  load  deflection  giving. 


_ (l  - g + u ) 

^ (1  -l  + a) 


1 - ^ ( 


^1-0  + U- 


s®  + 1 


^1  + |-  + u. , 


s®  + 1 


where,  in  this  expression  the  Laplace  variable  has  been  nondemens ionalized 
by  the  sling  load  frequency.  Placing  the  AEROCRANE  values  of  buoyancy  ratio 
and  sling  load  ratio  in  this  expression  yields 


0.167 


Z* 


0.67  s 


•5  2 


S®  + 1 


+ 1 
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The  coefficient  in  front  of  the  brackets  is  the  low  frequency  amplitude 

Z* 

found  above.  If  the  frequency  response  is  calculated  for  — = 0 and 

Zl 

O.U72  (the  value  used  in  the  analog  simulation),  the  resulting  character- 
istic is  shown  in  Figure  6 . At  low  frequencies  u^  is  in  phase  with 
(-9)  indicating  the  usual  result  that  a nose  down  pitch  gives  a forward 
translation.  The  corresponding  acceleration  is  small  as  noted  above. 

If  faster  inputs  are  applied,  specifically  at  frequency  ratios  higher 
than  — ' ^ then  the  phase  shifts  l80°  as  indicated  in  the  figure. 

That  is  for  rapid  control  motions,  in  fact  due  to  the  sling  load  dynamics, 
the  vehicle  translational  acceleration  response  to  attitude  changes  sign 
indicating  that  there  will  be  difficulties  in  rapidly  controlling  the 
vehicle  with  cyclic  pitch  through  attitude,  i.e.,  the  effect  may  appear 
to  the  pilot  somewhat  like  a control  reversal  when  the  pilot  is  primarily 
concentrating  on  controlling  position  over  the  ground.  This  effect  weald 
tend  to  be  more  significant  in  buoyantly  supported  vehicles  owing  to  the 
low  static  gain  of  the  acceleration  - attitude  transfer  function,  which 
would  tend  to  make  the  pilot  apply  more  rapid  control  motions,  is  discussed 
further  in  the  section  on  analog  simulation. 

Attaching  the  load  at  the  center  of  gravity  alleviates  this  effect 
slightly  since  it  increases  the  frequency  at  wtiich  phase  shift  occurs. 

In  fact,  this  appears  to  be  the  major  effect  of  sling  load  attachment 
point  on  the  vehicle  control.  The  effect  of  the  pitching  moment  produced 
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as  a result  of  the  location  of  the  sling  load  attachment  point  is  small 
compared  to  the  speed  stability  effect  as  noted  in  the  analog  simulation 
section  where  the  response  is  discussed  in  more  detail. 

Figure  7 shows  the  extent  of  the  frequency  band  where  phase  shift 
occ'ors  as  a function  of  sling  load  frequency  and  sling  load  attachment 
point.  Shortening  the  sling  load  cable  will  of  course  increase  the 
frequency  at  which  the  l80°  phase  shift  occurs. 

The  corresponding  transfer  function  for  a helicopter  assuming 
M,  =0.75  is  equal  to 


s _ 
gS 


(1  + 0.43  s^  + 1 


0.571  s^  + 1 


Without  buoyancy? the  low  frequency  response  amplitude  is  one.  The  sling 


load  attachment  point  location  has  less  effect  on  the  frequency  at  which 
phase  shift  occurs  as  shown  also  in  Fig’ure  7 . The  frequency  response 
of  the  helicopter  is  also  shown  on  Figure  6 . Owing  to  the  much  larger 
value  of  translational  acceleration  obtained  without  the  buoyant  force 
there  is  probably  less  tendency  to  make  rapid  control  motions  and  con- 
sequently less  tendency  to  operate  into  the  region  of  l80°  phase  shift. 

The  simple  model  employed  here  Indicates  that  this  low  value  of 
translational  acceleration  is  inherent  in  any  hovering  vehicle  which 
obtains  an  appreciable  portion  of  its  lift  from  buoyancy  and  that 
further  handling  qualities  studies  are  highly  desirable  to  fully  under- 
stand the  importance  of  this  marked  difference  from  a helicopter.  It 
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also  tends  to  indicate  that  short  sling  load  lengths  (high  sling  load 
frequencies)  are  desirable,  such  that,  the  pilot  will  not  encounter  the 
phase  shift  region. 

With  a translational  force  control  such  as  the  winglets,  this  effect 
is  not  present  since  the  control  directly  produces  translational  acceler- 
ation. 
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Figure  6 : Frequency  Response  Comparison  of  Translational 
Acceleration  to  Pitch  Attitude  Response  for 
AEROCRAIJE  and  Helicopter. 
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Figure  7:  Frequency  at  Which  Phase  Change  Occurs  in  Translational 
Acceleration  Pitch  Attitude  Response  as  a F’onction  of 
Sling  Load  Frequency  and  Attachment  Point . 
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NOK-DIMENSIONAL  EQUATIONS  OF  MOTION 


It  is  convenient  to  non-dimensionalize  the  equations  of  motion  of 
the  AEROCRANE  in  order  that  the  influence  of  vehicle  size  on  hovering 
stability  and  control  can  be  examined.  Further,  this  organization  of 
the  equations  of  motion  will  make  it  possible  to  readily  examine  the 
effects  of  various  parameters  of  interest  on  the  hovering  stability 
and  control. 

The  following  non-dimensionalization  is  employed.  Time  is  nondi- 
mensionalized  by  the  rotor  rotational  velocity  in  radians  per  second, 
n,  lengths  by  centerbody  radius,  R,,  and  forces  by  the  vehicle  weight 
W^.  The  following  notation  is  used  to  represent  nondimensional  quantities: 

( )^  = differentiation  with  respect  to  dimensionless 

time  (Qt). 

( ) = length  divided  by  centerbody  radius 

ou  = nondimensional  frequency,  ^ 

The  various  force  ratios  involved  in  the  equations  of  motion  are 
defined  as  follows : 

3 = buoyancy  ratio,  buoyant  force  divided  by  vehicle 

Fa 

weight,  — 
o 

u = sling  load  weight  ratio,  sling  load  weight  divided 

u 

by  vehicle  weight,  - — 

^o 

The  aerodynamic  thrust  to  weight  ratio  can  be  expressed  from  vertical 


i+3 


equilibrium  considerations  in  terms  of  the  ratios  p,  and  0 as 


T 


o 


= (1  + u - 0 ) 


The  apparent  mass  of  the  spherical  centerbody  is  equal  to  one  half 

(6) 

of  its  displaced  mass  and  consequently  is  equal  to 


m 2 
o 


It  may  be  noted  that  in  general,  this  relationship  between  the 
buoyant  force  ratio  and  the  apparent  mass  involves  the  assumption  that 
the  centerbody  is  filled  with  lifting  gas  (the  vehicle  is  not  equipped  with 
ballonets)  and  that  the  mass  of  the  lifting  gas  is  included  in  the  '/ehicle 
mass. 


With  these  definitions  it  is  possible  to  express  the  equations  of 
motion  for  the  hovering  dynamics  of  the  AEROCRANE  in  terms  of  the 

following  dimensionless  parameters  in  addition  to  0 and  p,  defined  above: 
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R,  = 
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ratio  of  centerbody  radius  to  rotor  radius , 

Froude  number  based  on  tip  speed  and  rotor  radius , 

W 


weight  coefficient. 


pnR^  (QR)' 


aerodynamic  moment  factor. 


2TtR' 


R 

g 


ratio  of  polar  moment  of  inertia  to  pitch/roll  moment  of  inertia 

W^Rs 

reference  frequency,  


kh 


m 2 _ 


' 3 L 


= uncoupled  sling  load  frequency,  — ^ 


Zl  n‘ 


Z*  = dimensionless  location  of  sling  load  attachment 
Za 

point,  — 

r = dimensionless  location  of  center  of  gravity 
o o ^ 

with  respect  to  center  of  buoyancy 

In  addition  to  these  parameters,  the  aerodynamic  stability  derivatives 
produced  by  the  rotor  enter  the  equations  of  motion.  In  non-dimensional 
form,  three  of  the  four  moment  derivatives  are  equal  (Appendix  B).  That 

2C„ 

or*  A 
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a CT  a a 
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a CT 


2C  „ 2C , 

U “-M. 

Expressions  for  the  two  derivatives  ■ and  ^ ^ are  given  in  Appendix  B. 
The  moments  are  nondimensionalized  by 


C M 5 L ~ 


M,L 


ottR®  (QR)®  R 


The  velocity  and  pitch  rate  are  nondimensionalized  by 


u q 

“ ' nS’  ■ n 


These  non-dimensional,  derivatives  are  a function  of  aerodynamic  thrust 
coefficient,  blade  lift  curve  slope,  rotor  solidity  and  centerbody  radius 
to  blade  radius  ratio.  Their  dependence  on  thrust  coefficient  is  given  in 
Figure  B-1  for  a solidity  typical  of  aEROCRAKE  designs.  The  theory  is 
de’/eloped  in  Reference  1 and  described  in  Appendix  B. 
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The  radius  of  gyration  of  the  vehicle  in  pitch,  k^,  can  be  introduced 
to  further  reduce  the  number  of  dimensionless  parameters . 


Consequently  two  of  the  parameters  given  above,  A and  can  be  expressed 


as 


A = 21!^  = 3 _g 

jf  ^ ..s  3 


W R, 


2 O 
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F,  k/x 


Alternate  forms  of  the  weight  coefficient  and  the  sling  load  frequency  are 
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' w 0 


pTTR2(nR)‘ 
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A 2 g 1 

^31=—^  = 


F,xz, 

Thus,  the  hovering  dynamic  stability  characteristics,  i.e.,  the  dimension- 
less natural  frequency  of  the  motion,  oun  , and  the  damping  ratio,  Q,  of  the 
AEROCEANE  with  sling  load  can  be  expressed  functionally  as 


•^n>C~F(F(,,0,j^  , k , x>  ^ y ^qJ  Z^,  Zi_5  ti) 

y ^ 


The  first  seven  parameters  are  associated  with  the  basic  design  of  the 
AEROCRAIIE, and  the  last  three  are  associated  with  the  sling  load  and  its 
attachment  point. 
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Consider  the  influence  of  size  on  the  dynamics  of  a series  of 
geometrically  similar  AEROCRAKES,  i.e.,  vehicles  with  the  same  values 
of  Xj  cr,  r^,  Z*,  Zf_.  To  examine  the  changes  in  the  dynamics  which  will 
occur  with  full  size  AEROCRAIffiS  in  contrast  to  the  model  of  Reference  1, 
four  parameters  must  be  considered:  the  buoyancy  ratio,  3;  the  radii  of 
gyration  and  ic^;  and  the  Froude  number  . 

The  buoyancy  ratio  is  largely  a design  choice.  Present  estimates 
of  the  size  of  this  ratio  for  a large  AEROCRANE  are  of  the  order  of  1.4. 
This  is  a somewhat  higher  value  than  the  dynamic  model  characteristics 
of  Reference  1 where  l.l4  < 3 < 1.19.  This  variation  would  be  expected 
as  certain  of  the  components  employed  in  a model,  e.g.,  the  power  plants 
tend  to  have  a higher  specific  weight  then  the  corresponding  ’units  on  a 
full  size  vehicle. 

Converting  the  thrust  to  weight  ratio  to  a thrust  to  sling  load 
(or  payload)  ratio  from  the  equation  above, 

-I_  = : + 1 - 9 

Wj  |_  |i 

Selecting  the  buoyancy  ratio  as  1.4-2  and  the  sling  load  ratio  as  0.84 
gives  a thrust  to  payload  ratio  of  0.5  as  is  characteristic  of  proposed 
AEROCRANE  designs. 

The  radii  of  gyration  might  be  expected  to  vary  from  the  dynamic 
model  characteristic  of  Reference  1 owing  to  different  construction 
techniques  which  would  be  employed  on  a full  size  AEROCRANE.  Largely 
the  radii  of  gyration  would  depend  upon  the  relative  mass  and  mass 
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distribution  in  the  blades  and  propulsion  units  compared  to  the  center- 
body  mass  distribution.  For  the  model  of  Reference  1,  the  complete 
moments  of  inertia  were  determined  to  be 
= 653  slug  ft^ 

l'  = 521  slug  ft^ 

includes  apparent  inertia  terms  from  the  blades  and  center  of  gravity/ 
center  of  buoyancy  spacing  amounting  to  71  slug-ft^.  Subtracting  this 
from  the  above,  = 450  slug  ft^.  About  60  percent  of  this  inertia 
arises  from  the  nacelles  and  blades  and  about  25  percent  from  the  instru- 
mentation oackage.  About  85  percent  of  the  nolar  moment  of  inertia  I 

* z 

is  contributed  by  the  blades  and  nacelles. 

The  radii  of  gyration  for  the  model  are 

iCy  =1.00  (without  apparent  inertia) 

k * =1.08  (including  apparent  inertia) 

y 

k = 1.21 
z 

It  would  be  expected  that  a full  size  AEROCRANE  would  have  similar 
inertial  characteristics  for  the  following  reasons.  A large  contri- 
bution to  the  pitch/roU  inertia  would  arise  from  the  gondola  crew 
station  which  would  take  the  place  of  the  instrumentation  package  on 
the  model.  The  relative  weight  of  the  blades  and  nacelles  would  tend 
to  be  less  than  the  model,  however,  the  propellers  would  represent 
relatively  more  contribution  to  the  inertia  since  the  wooden  propellers 


employed  on  the  model  were  very  light.  On  balance  it  might  be  expected 
that  the  dimensionless  radii  of  gyration  of  the  model  can  be  considered 
as  a first  approximation  to  the  characteristics  of  a full  size  AEROCRANE 
in  the  absence  of  a detailed  design.  For  reference,  Table  1 shows  the 
various  contributions  to  the  inertia  of  the  model  of  Reference  1.  The 
calculated  value  of  the  pitch  inertia  was  verified  by  experiment. 

Thus,  the  final  parameter  to  be  considered  is  the  Froude  number 

V = oh.  = 

" g Rg 

It  may  be  noted  that  this  parameter  is  a direct  measure  of  the  number 
of  g's  experienced  at  the  tip  of  a rotor  blade,  or  in  other  words,  the 
g field  in  which  the  power  plant  must  operate.  Taking  the  radius  as  a 
measure  of  the  size  of  the  vehicle,  the  angular  velocity  is  the  remaining 
parameter.  If  the  tip  speed  varies  with  size,  i.e.,  proportional  to  the 
radius,  what  may  be  considered  as  a series  of  Froude  scaled  vehicles 
results , and  the  rotor  RPM  would  vary  inversely  as  the  square  root  of 
the  linear  size  of  the  vehicle.  F'jrther,  the  dimensionless  dynamic 
characteristics  (natural  frequency  ratioed  to  rotor  angular  velocity  and 
damping  ratio)  would  be  invariant  with  size.  Therefore,  increasing  the 
size  of  the  AEROCRANE  in  this  fashion  would  only  result  in  a slower 
response  time  with  no  change  in  the  stability.  It  is  quite  reasonable 
to  expect  this  trend  with  size,  based  on  the  fact  that  the  rotor  thrust 
coefficient  with  full  payload,  would  be  invariant  with  size.  From  above 
Ct 

p-  = (1  + 0 - H) 
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or 


^ T 1,  y 3 I 


,(1^, 


Thus,  if  the  rotor  and  centerbody  geometry  (ct,  x)  invariant  with 

size  and  the  buoyancy  ratio  and  payload  ratio  have  been  selected  as 

(5)  ^ ^ 

noted,  the  average  rotor  blade  lift  coefficient  ^ (proportional  to  — ) 
is  determined  by  F,  . 

The  limitation  which  could  be  encountered  is  the  fact  that  the  tip 
speed  increases  with  size  to  maintain  the  Froude  number  constant.  However, 
this  vehicle  is  designed  with  a low  tip  speed  so  that  even  for  1 j,e  in- 
creases in  size  from  the  model  of  Reference  1,  only  moderate  tip  speeds 
result.  The  tip  speed  of  the  model  of  Reference  1 was  6l.7  feet  per 
second  nominally.  For  a large  AEROCRAHE  ten  times  the  size  of  this 
model  (rotor  radius  196.4  feet)  the  tip  speed  to  maintain  a constant 
Froude  number  and  consequently  the  same  blade  lift  coefficient  is  equal 
to  195  feet  per  second,  a very  moderate  value  compared  to  modern  heli- 
copter rotors . 

Thus,  it  is  considered  that  the  dimensionless  parameters  of  the 
model  described  in  Reference  1 are  reasonably  representative  of  a full 
size  AEROCRANE.  Table  II  lists  these  dimensionless  parameters  along 
with  those  obtained  in  a recent  design  study  for  a full  size  AEROCRANE 
with  a l6  ton  payload.  The  dimensionless  parameters  of  the  16  ton 
vehicle  are  used  for  the  study  of  the  influence  of  the  sling  load 
geometry  on  the  hovering  dynamics  discussed  in  another  section  of  this 
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report.  The  influence  of  the  center  of  gravity  position,  r^,  the  sling 
load  attachment  point,  Z^,  and  the  sling  load  length  are  examined. 

The  other  parameters  are  considered  fixed  at  the  values  gi’/en  in  Table 
II.  The  influence  of  size  of  the  vehicle  on  the  time  scale  of  the 
motion  is  obtained  by  determining  the  rotor  angular  velocity  from  the 
Froude  number  given  in  Table  II  and  the  radius  of  the  rotor  for  the  size 
AEROCRAWE  of  interest. 

If  an  AEROCRANE  equipped  with  winglets  is  considered,  then  an 
additional  parameter  is  introduced.  This  is  most  conveniently  expressed 
as  the  ratio  of  winglet  area  to  the  product  of  blade  chord  and  radius 
(cR).  If  it  is  assumed  that  the  chord  of  the  winglet  is  equal  to  the 
chord  of  the  blade  then  this  ratio  becomes  the  ratio  of  winglet  span  to 
blade  radius,  i.e.,  the  parameter  b^  introduced  in  Appendix  A. 
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TABLE  II 


^Parameters  shown  for  = 177  lbs.  only. 

(2 ) 

Umbilical  cable  weight  included  as  part  of  sling  load. 
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HOVERING  DYMMIC  STABILITY  CHARACTERISTICS 


OF  THE  AEROCRAJJE 

In  this  section,  the  dj-namic  stability  characteristics  of  the 
AEROCRAirE  in  hovering  are  described.  The  characteristic  roots  of  the 
AEROCRAtlE  are  discussed  as  a function  of  the  vehicle  center  of  gra'/ity 
location  and  sling  load  geometry.  Results  are  presented  based  on  the 
dimensionless  parameters  given  in  Table  II  for  a l6  ton  vehicle.  Thus, 
the  characteristic  roots  are  presented  divided  by  rotor  ang’ular  velocity, 
Q.  For  various  size  AEROCRANES  with  the  non-dimensional  characteristics 
given  in  Table  II,  the  influence  of  size  can  be  evaluated  by  determining 
the  angijlar  velocity  corresponding  to  the  gross  'veight  of  the  size  vehicle 
of  interest. 

The  formulation  of  the  dynamic  model  of  the  AERCCRAiJE  is  described  in 
detail  in  Appendixes  3 and  C.  A linearized  model  with  six-degrees -of - 
freedom  (vehicle  roll  and  pitch,  horizontal  and  lateral  velocity,  sling 
load  pitch  and  roll)  is  employed.  It  has  been  shown  in  Reference  1 that 
this  analytical  model  accurately  predicts  the  meas’xred  dynamic  stability 
and  response  characteristics  of  a Froude  scaled  dynamic  model,  described 
by  the  dimensionless  parameters  given  in  Table  II. 

Figure  8,  gives  the  characteristic  roots  and  eigenvectors  for  the 
reference  configuration.  This  reference  configuration  has  the  center 
of  gravity  located  below  the  rotor  hub  at  one-tenth  the  centerbody 
radius  (r^  = O.l),  the  sling  load  is  attached  one  centerbody  radius  below 
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the  center  of  gravity  (Z^  = l.O)  and  the  nondimensional  sling  load 
frequency  = 0.195  corresponding  to  a sling  load  length  of  125  feet 

for  the  proposed  l6  ton  design.  These  geometric  characteristics  were 
employed  for  the  analog  simixlation  studies  described  in  another  section 
of  this  report. 

There  are  five  oscillatory  modes  involving  significant  components 
of  all  the  variables  of  motion;  one  high  frequency  mode  which  is  primarily 
an  angular  motion  with  negligible  contributions  from  the  translational 
velocities;  two  intermediate  frequencies  corresponding  essentially  to 
the  sling  load  dynamics ; and  two  low  frequency  vehicle  modes . 

All  of  the  modes  of  motion  are  stable  for  the  reference  configuration. 
It  can  be  noted  from  the  eigenvector  configurations  that  each  lateral  vari- 
able (e.g.,  vehicle  roll  angle  and  vehicle  pitch  angle)  is  90  degrees  out 
of  phase  and  of  equal  amplitude  indicating  that  the  transient  modes  are 
circling  motions,  and  can  be  identified  by  the  direction  of  the  circling 
compared  to  the  direction  of  rotation  of  the  vehicle. 

The  regressing  vehicle  and  sling  load  modes  are  lightly  damped 
compared  to  the  advancing  modes.  The  rapid  angular  mode  and  the  advancing 
modes  are  well  damped. 

In  the  two  low  frequency  modes , the  sling  load  amplitudes  relative  to 
the  vehicle  are  approximately  eight-tenths  of  the  angular  motion  of  the 
vehicle,  indicating  that  the  sling  load  remains  relatively  near  the 
vertical.  In  the  rapid  angular  mode,  the  sling  load  amplitude  is 
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approximately  one  and  one  half  times  the  vehicle  angijlar  motion,  indi- 
cating excitation  of  the  sling  load  in  rapid  maneuvering.  In  the  sling 
load  modes  the  angular  motion  of  the  vehicle  is  approximately  one  quarter 
of  the  sling  load  motion  indicating  relatively  weak  coupling  to  the 
vehicle  motion. 

The  influence  of  vehicle  geometry  variations  on  the  characteristic 
roots  are  shox-m  in  Figures  9 through  10.  In  these  complex  plane  graphs 
the  modes  of  motion  are  shown  only  in  the  upper  or  Icurer  half  plane 
depending  upon  whether  they  are  regressing  (upper  half  plane)  or  advancing 
(lower  half  plane)  modes  for  clarity. 

Figure  9 shows  the  influence  of  vehicle  center  of  gravity  position. 
The  sling  load  attachment  point  is  maintained  at  a distance  of  one  center- 
body  radius  below  the  reference  center  of  gravity  position.  It  can  be 
seen  that  lowering  the  center  of  gravity  has  a stabilizing  effect  on  the 
low  frequency  modes  as  earlier  studies  indicated.  The  center  of  gravity 
position  has  a negligible  influence  on  the  sling  load  modes  and  a small 
effect  on  the  rapid  angular  mode. 

Figure  10  shows  the  influence  of  sling  load  attachment  point  (Z;^) 
and  sling  load  length  (or  sling  load  natural  frequency)  on  the  dynamics. 
Sling  load  natural  frequency  has  only  a minor  influence  on  the  low  fre- 
quency modes  and  the  rapid  angular  modes, but  of  course,  causes  marked 
changes  in  the  sling  load  modes  as  would  be  expected.  With  the  sling 
load  attached  at  one  centerbody  radius  below  the  center  of  gravity 
(Z^  = 1.0)  the  modes  of  motion  are  stable  for  all  sling  load  frequencies 
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investigated.  Attaching  the  sling  load  at  the  vehicle  center  of  gravity 
(Z;^  = O)  produces  a destabilizing  effect.  The  regressing  low  frequency 
mode  is  unstable.  This  is  a result  of  a loss  in  effective  pendulosity 
arising  from  placing  the  sling  load  attachment  point  at  the  center  of 
gravity  of  the  vehicle.  In  addition,  as  might  be  expected  the  regressing 
sling  load  mode  is  also  mildly  unstable.  There  is  only  a small  influence 
on  the  rapid  angular  mode  in  either  case. 

The  instability  which  arises  is  very  mild.  These  analytical  results 
assume  no  sling  load  damping.  It  would  be  expected  that  if  a small  amount 
of  sling  load  damping  were  included  these  modes  would  be  stable.  Thus, 
from  a dynamic  stability  viewpoint  no  advantage  is  to  be  gained  from  locating 
the  sling  load  attachment  point  at  the  vehicle  center  of  gravity.  This  is 
of  course,  a difficult  configuration  to  achieve  owing  to  the  vehicle  geometry. 

All  of  the  above  results  were  obtained  using  a digital  computer  program. 

The  real  time  simulation  studies  described  elsewhere  were  conducted  on  an 
analog  computer  using  the  vehicle  characteristics  of  the  reference  config- 
'uration  described  above.  In  addition,  the  analog  model  incorporated  terms 
to  account  for  the  nonlinear  effects  of  centerbody  drag  and  Magnus  force  as 
well  as  a small  amount  of  sling  load  damping  relative  to  the  vehicle  =0.1}. 
Winglets  were  incorporated  in  the  analog  model  as  well,  using  the  theoretical 
model  for  winglet  contributions  described  in  Appendix  A.  Figure  11  shows 
the  transient  response  of  the  vehicle  to  cyclic  pitch  and  winglet  deflection 
as  obtained  from  the  analog  computer.  The  cyclic  pitch  response  character- 
istics agree  with  the  digital  computer  solutions  for  the  characteristic  roots 
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described  above.  The  responses  are  shown  hovrever  in  real  time  (0  = 1.15 
rad/sec)  for  the  l6  ton  proposed  design.  Note  that  the  sling  load  motion 
in  response  to  a cyclic  pulse  is  largely  a result  of  the  angular  response 
of  the  vehicle. 

The  response  to  a winglet  input  is  shown  in  Fig’ure  lib.  Here  the 
time  history  exhibits  considerably  more  damping  as  a result  of  the  winglet 
contributions  to  the  stability  derivatives  of  the  vehicle.  In  particular 
the  vringlet  size  chosen  for  the  investigation  (bj^  = 0.25),  based  on  the 
considerations  described  in  Appendix  A result  in  a very  large  drag  damping, 
(X,^  = - O.Uoy)  and  a small  increase  in  speed  stability  (AM^  = O.COl).  The 
large  results  in  a well  damped  low  frequency  mode  as  may  be  seen  by 
comparing  Figure  lib  (with  winglets)  to  Figure  Ha  (without  winglets).  Two 
other  features  of  the  transient  response  of  the  AEROCRANE  to  a winglet 
input  may  be  noted.  The  rapid  winglet  pulse  input  results  in  considerable 
excitation  of  the  sling  load  as  a result  of  the  initial  translational 
acceleration  produced  by  the  winglet.  Since  the  sling  load  has  a signi- 
ficant mass  compared  to  the  mass  of  the  vehicle  this  results  in  a coupling 
into  the  horizontal  velocity  as  may  be  noted  from  the  time  history.  Also 
as  the  vehicle  translates  fonvard  in  response  to  the  winglet  input,  a 
significant  pitch  up  occurs  as  a result  of  the  speed  stability  of  the 
rotor  resulting  in  a considerably  lower  steady-state  velocity  from  the 
winglet  input  than  would  be  expected  if  the  vehicle  remained  level.  This 
result  indicates  that  it  would  be  desirable  to  incorporate  attitude 
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feedback  to  cyclic  to  minimize  the  coupling  arising  from  translation 
if  the  winglet  controls  are  employed. 

Gusts  inputs  were  incorporated  in  the  analog  simulation  as  well. 
Unfortunately,  no  experimental  data  are  available  on  the  gust  response 
of  the  AEROCRATfE.  In  the  absence  of  data,  a simplified  model  was  em- 
ployed to  estimate  the  influence  of  wind  gusts  on  vehicle  motion.  Gust 
inputs  were  treated  as  equivalent  to  a change  in  the  translational 
velocity  of  the  vehicle  and  no  attempt  was  made  to  account  for  the  in- 
fluence of  the  spatial  velocity  distribution  of  a gust  field  on  the  re- 
sponse of  the  vehicle. 

Proper  modelling  of  the  gust  response  of  aEROCRANE  vehicles  is  a 
very  significant  aspect  of  precision  hovering  and  experimental  data  are 
essential  for  future  simulation  studies. 

Figure  12  shows  the  response  of  the  l6  ton  payload  configuration  with 
and  without  winglets  to  a step  gust  of  5 fps  and  a random  gust  with  an 
RMS  value  of  5 fps. 

A step  gust  input,  while  not  physically  realistic,  does  provide  a 
good  basis  for  comparison  of  the  gust  response  of  the  vehicle  with  and 
without  winglets. 

It  can  be  seen  by  reference  to  Figures  12a  and  12b  which  show  the 
response  to  a step  gust,  that  with  winglets,  the  initial  translational 
velocity  response  is  considerably  faster.  There  is  an  initial  horizontal 
acceleration  due  to  the  presence  of  the  winglets  and  the  vehicle  reaches 
the  steady-state  velocity  of  5 fps  in  a somewhat  shorter  time.  With 
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winglets,  the  long  period  dynamics  are  well  damped  while  in  the  winglet  off 
case,  the  lightly  damped  oscillation  persists.  Figures  12c  and  12d  show  the 
uncontrolled  response  of  the  vehicle  with  and  without  winglets  for  a random 
gust  disturbance  with  a 5 fps  RMS  level.  The  amplitude  of  the  response  of 
the  vehicle  is  similar  in  both  cases  however,  with  winglets  there  is  a notice- 
able increase  in  the  high  frequency  content  of  the  response  especially  in 
the  horizontal  velocity  time  history.  The  presence  of  the  winglets  results 
in  forces  acting  on  the  vehicle  which  are  linearly  proportional  to  velocity 
and  in  addition,  the  large  value  of  noted  above  results  in  a comparatively 
fast  time  constant  in  the  translational  response  of  the  vehicle  which  is  not 
present  when  the  winglets  are  absent.  This  increase  in  the  high  frequency 
response  to  gusts  would  indicate  that  the  vehicle  equipped  with  winglets 
would  experience  higher  structural  loads  in  turbulent  air. 

The  next  section  of  this  report  describes  the  results  of  some  real  time 
simulation  studies  of  the  AEROCRANE  with  the  vehicle  dynamic  characteristics 
illustrated  in  Figures  11  and  12. 
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Figure  8a:  Eigenvectors  for  Reference  Case.  Low- 
Frequency  Modes. 
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Figure  8b:  Eigenvectors  for  Reference  Case.  Sling  Load 
and  High  Frequency  Modes. 
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Figure  9a : Influence  of  Center  of  Gravity  Location . 
Low  Frequency  Modes, 
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Figure  9b:  Influence  of  Center  of  .Gravity  Location.  Sling  Load 
and  High  Frequency  Modes. 
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Figure  lOst  Influence  of  Uncoupled  Sling  Load  Frequency  and  Sling 
Load  Attachment  Point  on  Low  Frequency  Modes. 
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Figure  10b:  Influence  of  Uncoupled  Sling  Load  Frequency  and  Sling  Load 
Attachment  Point  on  Sling  Load  and  High  Frequency  Modes. 
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Figiire  11a:  Transient  Response  of  l6  Ton  AEROCRANE  to  Cj'clic 
Pulse  Input. 
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Figure  12a:  Gust  Response  of 
5 fps  Step  Gust. 
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Ton  AEROCRANE.  Winglets  Off 


SIMULATION  STUDIES 


An  investigation  of  the  handling  qualities  of  the  AEROCRANE  in 
hovering  flight  was  conducted  using  an  analog  computer  to  model  the 
vehicle  dynamics  in  and  near  hovering  flight.  The  equations  of  motion 
programmed  are  described  in  the  Appendixes.  A six-degree -of -freedom 
model  was  employed  including  four-degrees-of-freedom  for  the  basic 
vehicle  and  two-degrees-of-freedom  for  the  sling  load.  A random  gust 
input  was  employed  to  disturb  the  vehicle  for  certain  of  the  experiments. 

Only  the  reference  configuration  was  examined,  and  the  vehicle  parameters 
employed  were  those  of  a l6  ton  proposed  AEROCRANE  design  given  in  Table  III. 

The  basic  objective  of  this  investigation  was  to  obtain  a compari- 
son of  the  handling  of  the  AEROCRANE  in  hover  with  the  heavy  lift  heli- 
copter (HLH)  studies  reported  in  Reference  3,  and  in  addition,  to  compare 
the  control  of  the  AEROCRANE  with  cyclic  pitch  to  control  with  winglets. 

Consequently,  two  tasks  were  examined;  hovering  over  a spot  with 
the  operator  attempting  to  maintain  a position  in  space  in  the  presence 
of  random  gusts  with  an  RMS  level  of  5 fps;  and  precision  maneuvering, 
that  is  flying  from  one  location  to  another. 

It  should  be  noted  that  the  set-up  involved  was  a very  simple  one. 

The  longitudinal  and  lateral  displacements  of  the  vehicle  in  space  were 
fed  into  an  x-y  plotter  and  the  operator  attempted  to  perform  the  task 
with  only  information  on  the  position  of  the  vehicle  in  space  through  a 
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proportional  control  stick.  No  simulated  cockpit  was  employed  and  the 
operator  did  not  have  the  benefit  of  any  cues  other  than  the  motion  of 
the  vehicle  as  perceived  cn  the  x-y  plotter.  It  is  felt  that  this  par- 
ticular experimental  set-up  tends  to  make  cyclic  control  in  hovering 
inherently  difficult. 

In  general,  with  cyclic  control  it  was  necessary  to  employ  attitude 
feedback  to  be  able  to  achieve  any  reasonable  level  of  precision  control 
as  would  be  expected  given  the  task.  The  attitude  feedback  employed  was 
of  a form  found  to  be  desirable  as  a consequence  of  the  experimental 
results  of  Reference  1.  That  is,  feedback  of  the  form 
Ai3  = -K|^(9  -0) 

3is  = K*(e  + 0) 

were  employed  with  a gain  0.3  . Basically  it  was  found  that 

using  winglet  control  either  of  the  tasks  was  relatively  easy  to  perform, 
while  with  cyclic  control  both  tasks  were  considerably  more  difficult. 

The  addition  of  attitude  feedback  appeared  to  have  only  a small  effect 
on  the  control  task  with  winglet  control.  However,  neither  of  the 
tasks  could  be  accomplished  with  cyclic  control  without  attitude  feed- 
back. 

First  consider  the  task  of  hovering  over  a spot.  Figure  13  shows 
a comparison  of  the  time  history  of  hovering  over  a spot  with  winglet 
control  to  that  with  cyclic  control.  The  gust  level  employed  in  these 
two  oases  was  5 fps  RMS.  The  operator  was  able  to  keep  the  vehicle. 


much  of  the  time,  x^^ithin  a one  foot  square  using  vinglet  control  and  only 
ocassionally  was  a two  foot  square  exceeded.  The  total  time  of  the  exper 
ment  shown  was  five  minutes  and  no  attitude  feedback  was  employed  (recall 
that  the  vehicle  has  an  inherent  attitude  feedback  due  to  the  center  of 
gravity/center  of  buoyancy  spacing).  Thus,  using  winglet  control  it 
appears  that  the  basic  vehicle  with  no  automatic  stabilization  can  be 
very  effectively  hovered  over  a spot.  In  fact,  indications  from  this 
experiment  indicate  the  the  winglets  give  precision  control  of  a level 
similar  to  that  achieved  with  a sophisticated  automatic  control  system 
in  the  KLH  simulation  of  Reference  3 as  shown  in  Figure  l4. 

With  cyclic  control  the  situation  was  quite  different.  It  was  not 
possible  to  maintain  a position  within  the  scale  of  the  x-y  plotter 
(±  5 feet)  for  any  length  of  time  with  the  basic  vehicle.  With  attitude 
feedback  it  was  possible  to  maintain  position  within  a 10  foot  square 
for  three  minutes  as  shown  in  Figure  13.  This  appears  to  be  somewhat 
poorer  performance  than  indicated  for  the  HLH  with  attitude  hold  as 
shown  in  Figure  l4. 

Now  consider  the  other  task,  that  of  precision  maneuvering  from 
point  to  point  at  low  speeds.  The  operator  commenced  flying  at  the 
center  of  the  x-y  plotter  and  then  flew  to  a one  foot  square  box  at 
one  location  and  then  on  to  the  next  as  shown  in  Figure  15 . No  gust 
inputs  were  employed  in  this  task.  Again  with  winglet  control  the 
task  was  relatively  easy.  No  significant  difference  was  noted  between 
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the  cases  with  and  without  attitude  feedback.  A relatively  straight 
path  between  boxes  was  flown  and  there  was  no  difficulty  stopping 
and  remaining  within  the  box  as  shown  in  Figure  15 . This  behavior  is 
quite  similar  to  that  indicated  in  Reference  3 for  the  HLH  with  the 
high  gain  velocity  hold  system  as  shown  in  Figure  lU  . 

With  cyclic  control  again  the  task  proved  impossible  without 
attitude. feedback.  With  attitude  feedback  a sample  trace  is  shown 
in  Figure  15.  The  task  was  still  very  difficult  and  a considerable 
amount  of  practice  was  required  (approximately  20  minutes)  to  achieve 
the  result  indicated  in  the  figure.  It  appeared  to  the  operator  that 
precision  translation  must  be  accomplished  either  at  a very  low  speed 
or  a very  "high"  speed.  The  operator  noticed  significant  inputs  from 
the  sling  load  motion  which  complicated  the  control  task.  The  diffi- 
culty with  this  task  is  clearly  related  to  the  problem  area  discussed 
in  the  section  on  sling  load  dynamics  where,  when  low  frequency  control 
is  attempted  only  a small  translational  acceleration  is  achieved.  As 
higher  frequency  inputs  are  employed  the  translational  acceleration 
decreases  due  to  sling  load  motion  and  ultimately  the  phase  of  the 
attitude/acceleration  response  can  reverse.  It  should  be  noted  that 
this  characteristic  appears  to  be  rather  fundamental  to  precision 
control  of  vehicles  which  obtain  an  appreciable  amount  of  their  lift 
from  buoyancy  as  noted  in  another  section. 

In  order  to  ascertain  in  a simple  fashion  the  importance  of  the 
sling  load  dynamics,  the  sling  load  damping  ratio  was  increased  to  0.5. 
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The  operator  noted  that  the  sling  load  didn't  tend  to  "throw  the  vehicle 
around"  as  much,  and  it  was  considerably  easier  to  achieve  the  desired 
result  of  flying  from  place  to  place  as  shown  in  Figure  15- 

It  should  also  be  noted  that  in  all  of  these  experiments  the 
operator  was  controlling  the  center  of  gravity  motion  of  the  vehicle. 

The  sling  load  motions  corresponding  to  two  typical  runs  are  shown  in 
Figure  15  indicating  that  there  is  considerable  sling  load  motion.  In 
addition,  another  feature  of  the  hovering  control  of  the  AEROCRANE  which 
may  be  of  considerable  importance  is  the  influence  of  the  location  of 
the  pilot  in  the  actual  vehicle  with  respect  to  the  center  of  gravity 
of  the  vehicle,  with  respect  to  the  hovering  control  problem.  Some 
discussion  of  considerations  in  this  regard  can  be  found  in  Reference 
10.  It  was  not  possible  with  the  simple  simulator  employed  in  these 
studies  to  investigate  this  point. 


77 


TABLE  III 


PHYSICAL  PARAMETERS  OF  l6  TON  AEROCRAJ^IE  USED  IN 
ANALOG  SIMULATION  STUDY 


1.  ) Overall  Vehicle 

= 38,000  lbs. 

Fb  = 54,000  lbs. 

l'  = 5.40  X 10®  slug  ft^ 

I^  = 6.75  X 10®  slug 

= 5.9  ft. 

2.  ) Rotor 


R = 

136.5  ft. 

c = 

15.25  ft. 

n = 

11  rpm  (1.15 

a = 

5.73 

b 

li 

T 

16,000  lbs. 

Centerbody 

Rb  = 

59  ft. 

Co  = 

0.78 

C u M ~ 

0.30 

Sling  Load 

= 

32,000  lbs. 

Za  = 

59  ft. 

Zu  = 

125  ft. 

II 

0.10  . 

Winglet* 

Cw  = 

15.25  ft. 

bw  = 

34.13  ft. 

aw  = 

4.6 

*Note  that  for  the  simulation  studies  involving  cyclic, 
the  winglets  are  removed  from  the  model. 
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Figure  13:  Comparison  of  Hovering  Hold  in  Cast  with 
Winglet  Control  and  Cyclic  Control. 

16  Ton  AEROCRAHE. 
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Figure  lU:  Simulation  Results  for  Heavy  Lift  Helicopter 
(Reference  3)« 
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Figure  15a:  Precision  Maneuvering  with  Winglet  Control. 
16  Ton  AEROCRANE. 
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COMPARISON  OF  AEROGRAMS  HOVERING  DYNAMICS  WITH 


COirVENTIONAL  HELICOPTERS 

It  is  instructive  to  compare  the  hovering  stability  and  control 
characteristics  of  the  AEROCRANE  with  conventional  helicopters  of  equiv- 
alent size.  Reference  5 presents  estimates  of  the  stability  derivatives 
of  various  large  proposed  tandem  rotor  helicopter  designs  with  gross 
weights  varying  from  25  tons  to  100  tons.  Reference  U presents  the 
hovering  stability  and  control  derivatives  for  a 20  ton  single  rotor 
helicopter.  Thus,  to  include  single  rotor  in  this  discussion,  an  es- 
timate of  ihe  variation  of  these  derivatives  with  size  of  the  vehicle 
must  be  made . As  indicated  in  Reference  U , the  variation  of  the 
stability  and  control  derivatives  of  a conventional  helicopter  with 
gross  weight,  will  be  somewhat  different  than  the  variations  indicated 
for  the  AEROCRANE  primarily  as  a result  of  the  fact  that  for  various 
technical  reasons  all  helicopters  operate  at  about  the  same  tip  speed 
and  consequently  as  the  vehicle  size  is  increased  the  tip  speed  is 
maintained  constant  rather  than  increasing  with  size  as  Froude  scaling 
would  indicate.  The  consequence  of  this  on  the  variation  of  the 
important  stability  derivatives  can  be  determined  in  the  following  way. 
Consider  the  variation  of  the  three  stability  derivatives  largely 
responsible  for  the  longitudinal  hovering  dynamic  characteristics  of  a 
vehicle;  the  pitch  damping,  M^,  the  speed  stability  M^,  and  the  control 
effectiveness  NL  . Expressed  in  non-dimensional  form: 


83 


AD-A054  2ai 


uNCLASSXFZCO 


PHXNCCTON  UNXV  N J OCPT  Of  AEROSPACE  AND  MECHANICAL— ETC  f/O  1/3 
A STUOr  Of  THE  PRECISION  lOVER  CAPABILITIES  Of  THE  AEROCRANE  HY— ETC(U) 
fCB  7S  H C CURTlSSf  H f PUTMAN*  R M MCKILLIP  N62269-77-C-0074 


AMS-Tfl-1949 


NAOC-7A341-30 


^ .pttr^  (QR)^  R,  1 


= (Cis!,Mlls,  ^ c_ 


y 


QR  m 
M> 


^ pnR^  (QR)^  R g 


5 C j and  C will  depend  only  upon  the  geometry  of  the  helicopter. 

IHa  ni  QL 

1 M.  0 

Thus,  if  a series  of  geometrically  similar  hovering  vehicles  are  considered 
the  variations  in  the  three  dimensional  derivatives  with  size  can  be  ex- 
pressed as 


R*^  Vt 


M 


M 


R®  Vt 


y 

R^  Vt^ 


where  V x = QR 

For  the  conventional  helicopter  it  may  be  assumed  that  V x is  independent 
of  size  (Reference  U),  consequently 


M ~ ^ 
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Defining  the  ratio  of  the  gross  weight  of  one  helicopter  to  another  of 
similar  geometry  as 


Reference  4 indicates  that  as  the  helicopter  gross  weight  is  increased, 
the  linear  dimensions  grow  as  the  square  root  of  the  weight  ratio 


The  moments  of  inertia  increase  with  gross  weight  with  a power  somewhat 
less  than  the  gross  weight  ratio  squared. 


Substituting  these  ratios  into  the  above  expressions  it  is  found  that 


~ w 


.1 


.4 

w 


The  pitch  damping  is  approximately  independent  of  size.  The  speed 
stability  and  the  pitch  control  power  decrease  with  increasing  gross  weight 
varying  approximately  as  the  square  root  of  the  gross  weight.  Table  I of 
Reference  4 indicates  that  these  simple  relationships  give  the  dominant 
trends  in  the  pitch  derivatives.  The  roll  derivatives  vary  somewhat 
differently  owing  to  the  smaller  increase  in  roll  inertia  with  weight. 


85 


For  the  AEROCRAKE  a somewhat  different  dependence  is  obtained.  In 
particular,  since  the  tip  speed  of  the  AEROCRANE  is  low,  the  variations 
in  these  quantities  are  somewhat  different.  It  has  been  assumed  that  a 


series  of  Froude  scaled  vehicles  is  being  considered,  such  that 


There  is  no  difficulty  encountered  with  increasing  the  tip  speed  as 

required  by  Froude  scaling,  and  in  this  case  the  linear  dimensions 

grow  as  the  cube  root  of  the  weight.  The  average  rotor  blade  lift 
Ct 

coefficient  (— ) is  independent  of  size.  Note  that  if  constant  tip 
speed  is  assumed,  the  rotor  radius  grows  as  the  square  root  of  the 
gross  weight  to  maintain  the  same  blade  lift  coefficient.  In  this 
case  the  following  derivative  variations  result. 
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Thus  in  Froude  scale  all  of  the  derivatives  are  reduced,  the  largest 
reduction  occurring  in  the  speed  stability.  It  may  be  recalled  that 
in  the  Froude  scaled  case  there  is  no  change  in  the  dimensionless 
vehicle  dynamics,  that  is,  only  the  time  scale  of  the  motion  changes. 

Table  IV  lists  the  derivatives  for  various  size  vehicles.  The 
values  for  the  tandem  helicopter  are  taken  from  Reference  U.  The 
values  for  the  single  rotor  are  taken  from  Reference  10 for  a 20  ton 
vehicle  (S-65)  and  then  are  scaled  to  a 100  ton  vehicle  using  the 
constant  tip  speed  scaling  described  above.  The  AEROCRANE  derivatives 
are  based  on  the  Froude  scaling  given  above. 

There  are  some  other  fundamental  differences  between  the  aEROCRAWE 
and  helicopters  to  be  noted  before  discussing  the  influence  of  these 
various  derivatives  on  the  motion.  First  it  may  be  noted  that  the 
influence  of  the  speed  stability  on  the  dynamics  will  be  considerably 
less  for  the  AEROCRANE  than  for  the  helicopter  owing  to  the  fact  that 
the  vehicle  is  supported  in  part  by  buoyancy.  The  influence  of  the 
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attached  to  the  vehicle  or  in  other  words,  that  9 = 9^  which  is  approxi- 
mately true  in  the  long  time  dynamics.  Consequently  this  factor,  which 
appears  multiplied  times  in  the  characteristic  equation,  is  about 
25  percent  of  g for  the  parameters  assumed  for  the  AEHOCRAilE  and  conse- 
quently there  is  less  influence  of  on  the  characteristic  roots . 

Another  way  of  understanding  this  effect  is  that  with  a buoyant  force, 
less  translational  acceleration  is  produced  by  the  AEROCRANE  as  a result 
of  pitch  and  consequently  there  is  a smaller  velocity  excursion  associated 
with  pitch  and  consequently  less  impact  of  on  the  dynamics. 

One  further  aspect  of  the  AEROCRANE  tends  to  improve  the  inherent 
dynamics  of  the  AEROCRANE  compared  to  hovering  helicopter  and  that  is 
the  fact  that  the  center  of  gravity  is  located  below  the  center  of 
buoyancy  and  consequently  there  is  some  inherent  attitude  stability  which 
is  no  present  on  a conventional  helicopter  without  attitude  feedback. 

Certain  other  considerations  also  result  in  differences  in  the 

dynamic  behavior.  Consider  the  short  time  steady-state  degree-of-freedom 

motion  with  a step  input  in  cyclic.  Using  the  dimensionless  aerodynamic 

derivatives  for  the  AEROCRANE  from  Appendix  B 

- C = C 
m-  m. 
a \s 

Consequently  when  the  vehicle  has  attained  a steady  pitch  rate 
M ^ + M,  Ai,  = 0 

‘1  IS 

Using  the  relationships  given  above  for  the  dimensional  derivatives 
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taken  with  the  equality  of  the  non  dimensional  derivatives,  the  equation 
above  becomes 


9 — + n A j 


that  is,  the  short  time  angular  rate  developed  by  a cyclic  input  is 
directly  proportional  to  the  rotor  angular  velocity.  Consider  now  a 
simple  model  for  the  single  rotor  helicopter.  Assuming  that  the  total 
moment  applied  to  the  single  rotor  helicopter  is  proportional  to  longi- 
tudinal flapping,  the  equation  for  the  flap  angle  can  be  written  as 


A steady  angular  rate  implies  that  there  is  no  moment  applied  to  the 
helicopter,  consequently  Aa^j  = 0.  Therefore 


^ = 


yQ 


>13 


Y is  the  order  of  10  for  the  single  rotor  helicopter  discussed  in 
Reference  10,  consequently 

9 = " § ^ Bis 


and  a result  of  similar  form  to  the  AEROCRANE  is  obtained.  This 
approximate  result  also  applies  to  the  single  rotor  helicopter  and 
the  tandem  helicopter  in  roll.  For  the  tandem  helicopter  in  pitch 
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a somewhat  different  expression  results  since  both  control,  which  is 
obtained  from  differential  collective  pitch,  and  damping  arise  from 
thrust  variations.  Assuming  a tandem  helicopter  with  the  center  of 
gravity  located  midway  between  the  two  rotors,  the  dimensionless 


derivatives  are  approximately  given  by 


C 

nu 

q 


0 

0 


i it  2 
n 3 


9 


= 0 
o 


Thus , it  can  be  seen  in  al 1 three  cases  that  the  fundamental  parameter 
determining  the  angular  steady-state  rate  for  a given  control  input  is 
the  rotor  RPM.  Consequently  the  AEROCRANE  with  its  lower  tip  speed  will 
tend  to  develop  a smaller  pitch  rate  for  a given  physical  control 
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deflection  than  the  conventional  helicopter,  i.e.,  a larger  cyclic 
deflection  with  stick  motion  will  tend  to  be  desirable. 

The  l6  ton  payload  AEROCRAHE  has  a rotor  RIM  of  11  (1.15  rad  per  sec) 
whereas  the  40  ton  gross  weight  tandem  has  a rotor  RPM  of  I9.8  rad  per  sec 
indicating  a considerably  larger  pitch  rate  for  a given  angular  deflection 
of  the  control. 

I 

The  actual  derivatives  give  a steady  rate  due  to  cyclic  of  7.0  degress 
per  sec/inch.  To  achieve  this  value  on  the  AEROCRAKE  would  require  for  a 1 

rotor  RPM  of  11,  a gearing  of  the  order  of  6 degrees  of  cyclic  per  inch  of 
stick.  The  implications  of  this  are  discussed  below  in  terms  of  handling 
qualities. 

Comparison  of  the  AEROCRANE  damping  and  control  sensitivity  in  hovering 
with  the  tandem  helicopter  characteristics  given  in  Reference  4 and  the 
single  rotor  data  given  in  Reference  10  is  shown  in  Figure  16  . The  single 
rotor  data  given  in  Reference  10  is  for  a 20  ton  vehicle  (S-65)  for  the 
roll  axis  only.  To  obtain  pitch  axis  characteristics  it  is  assumed  that 
the  ratio  of  pitch  inertia  to  roll  inertia  is  4.5  and  that  the  longitudinal 
stick  gearing  is  2.5  times  the  lateral  stick  gearing.  The  characteristics 
of  a 100  ton  single  rotor  were  obtained  by  scaling  the  20  ton  vehicle 
using  the  parameter  variations  with  gross  weight  given  in  Reference  4 and 
described  above.  Reference  10  indicates  a similar  variation  in  roll 
inertia  for  the  single  rotor  as  Reference  4 indicates  for  the  tandem. 

For  comparison  purposes,  the  control  gearing  of  the  AEROCRANE  was  taken 
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to  be  three  degrees  of  cyclic  per  inch  of  control  deflection.  This 
gearing  is  chosen  for  comparison  purposes  only. 

Also  shown  on  Figure  l6  are  the  3-|  pilot  rating  boundaries  deter- 
mined by  the  handling  qualities  experiments  reported  in  Reference  11  for 
a 15  ton  single  rotor  helicopter. 

It  should  be  kept  in  mind  in  the  following  discussion  that  this  is 
the  largest  helicopter  for  which  handling  qualities  data  are  available. 

The  Military  Specification  for  Helicopter  Flying  and  Ground  Handling 
Charactei istics  (MIL-H-85OIA)  indicates  a downward  trend  of  damping  and 
control  sensitivity  with  size  and  it  is  concluded  in  Reference  11  that 
the  results  presented  there  generally  support  the  desirable  trends  indicated 
in  KLL-H-d^ClA.  It  should  be  noted  that  the  criteria  of  MIL-H-85OIA  were 
derived  from  flight  experiments  on  a 2.5  ton  helicopter.  It  should  also 
be  noted  that  the  AEROCRAKE  differs  in  two  respects  from  the  conventional 
helicopters : 

1. )  Center  of  gravity/center  of  buoyancy  spacing  provides  an 

inherent  attitude  stability  and  L^)  which  will  in- 
fluence judgements  on  handling  qualities. 

2. )  The  translational  acceleration  of  an  AEROCRANE  per  degree 

of  tilt  is  about  one-sixth  of  that  of  a helicopter. 

The  presence  of  seme  attitude  stability  is  certainly  in  a favorable 
. iction.  It  is  difficult  to  assess  the  impact  of  the  translational 
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acceleration -attitude  relationship. 

Figure  l6  shows  then  the  comparison  of  these  vehicles.  Generally 
the  AEROCRAIIE  characteristics  place  it  on  the  left  of  the  graphs  of 
damping  vs.  control  sensitivity  compared  to  the  helicopters.  The  AERO- 
CRANE  in  general  has  a very  high  level  of  angular  damping  about  both 
axes  and  a low  control  sensitivity  for  the  assumed  stick  gearing.  It 
does  however  lie  quite  close  to  the  visual  flight  S'i  pilot  rating 
boundaries.  Its  characteristics  are  essentially  on  the  pilot  rating 
boundary  in  the  roll  case  and  somewhat  outside  the  pitch  boundary.  The 
characteristics  of  the  AEROCRANE  relative  to  the  helicopters  is  primarily 
a result  of  the  low  rotor  angular  velocity  described  above. 

Figure  17  shows  a comparison  of  the  AEROCRANE  damping  and  control 
sensitivity  characteristics  compared  to  the  boundary  given  in  Reference 
11  along  with  that  of  Reference  12  (MIL-H-85OIA) . It  can  be  seen  that 
both  AEROCRANE  vehicles  examined  essentially  meet  the  specifications  of 
MIL-H-85OIA.  The  relationships  given  in  MIL-H-85OIA  were  evaluated  for 
a 35  ton  gross  weight  vehicle  in  the  case  of  the  I6  ton  payload  AEROCRANE 
and  a 109  ton  gross  weight  vehicle  for  the  50  ton  payload  configuration. 
Again  it  should  be  emphasized  that  applying  the  criteria  of  Reference  12 
to  vehicles  of  this  gross  weight  represents  a considerable  extrapolation 
from  the  vehicle  on  which  the  data  are  based  (2.5  ton)  as  well  as  the 
vehicle  of  Reference  11  (15  ton).  As  an  example,  the  minimum  values  of 
the  angular  damping  derivatives  given  by  Reference  12  are 


- 18  I ' 

X 


.3 
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For  the  16  ton  and  50  ton  AEROCRANE  these  expressions  give; 


for  the  16  ton  payload  vehicle 
M = - 0.076  per  sec 

q. 


Lp  = - 0.172  per  sec 


and  for  the  50  ton  payload  vehicle 


Mg  = - 0.0^3  per  sec 


L = - 0.097  per  sec 
P 


These  required  levels  of  damping  indicated  in  Reference  12  are 
essentially  more  than  a factor  of  10  less  than  the  values  estimated 
for  the  AEROCRAJJIE.  As  can  be  seen  by  this  comparison  as  '■veil  as 
from  the  Figures,  the  AEROCRANE  possesses  a very  high  level  of 
angular  damping  and  so  it  appears  difficult  to  draw  detailed  con- 
clusions from  this  comparison  with  other  vehicles  as  well  as  with  the 
specifications.  Increasing  the  stick  gearing  above  the  value  assumed 
for  the  AEROCRANE  will  place  the  vehicle  well  within  the  criteria  of 
Reference  12,  however  there  is  a significant  extrapolation  involved 
both  as  a result  of  the  very  large  size  of  the  vehicle  and  the  ■tnvo 
unique  characteristics  noted  abo-ve  to  assume  that  meeting  the  speci- 
fications of  MIL-H-85OIA  will  produce  satisfactory  handling  qualities 
in  hovering. 

The  roll  and  pitch  axes  of  the  AEROCRANE  are  coupled  as  they  are 
on  a single  rotor.  The  tandem  essentially  has  little  or  no  pitch/roU 
coupling.  The  coupling  levels  of  the  AEROCRANE  are  similar  to  those 
indicated  in  Table  III  for  the  single  rotor.  In  both  cases,  the  coupling 
derivatives  (L^^,  M^)  are  of  a similar  size  to  the  damping  derivatives 
(«,>  i-p)- 
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The  speed  stability  and  dihedral  effect  of  the  AEROCRANE  are 
larger  than  the  values  for  the  tandem  indicating  a larger  gust  sensi- 
tivity. However,  recall  that  the  .AEROCRAlffi  will  have  less  translation 
as  a result  of  pitching  so  that  it  is  difficult  to  completely  compare 
the  implications  of  this  difference.  The  speed  stability  of  the  AERO- 
CRAKE  is  about  three  times  that  of  the  tandem,  however  pitch  of  the 
AEROCRAME  gives  about  one-sixth  the  translational  acceleration  of 
the  tandem  so  the  net  effect  of  a horizontal  gust  acting  through  the 
speed  stability  on  the  translational  acceleration  is  of  a similar  size. 

In  the  case  where  winglet  control  is  used  it  is  more  difficult  to 
find  criteria  regarding  handling  qualities.  Reference  13  indicates 
that  for  aircraft  which  use  direct  side  force  as  a primary  lateral 
positioning  device,  the  direct  side  force  control  should  be  sufficiently- 
powerful  to  obtain  lateral  acceleration  values  between  0.08  and  0.12  g 
in  wings-level  flight.  The  winglet  size  employed  in  the  analog  simulation 
studies,  had  a control  sensitivity  of  = l.l6  ft  per  second  squared 


per  degree,  indicating  that  a winglet  deflection  of  3 would  produce  a 
satisfactory  level  of  translational  acceleration.  Again,  it  should  be 
noted  that  this  value  was  derived  from  handling  qualities  studies  on  a 
2 ton  aircraft.  It  is  also  noted  in  Reference  13  that  "information  is 
needed  for  side-force  controls  on  attitude-stabilized  aircraft  before 
generalizations  to  other  aircraft  can  be  made". 

The  analog  simulation  studies  reported  here  used  a total  winglet 

Q 

deflection  of  3 resulting  in  a maximum  acceleration  equal  to  the  value 
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indicated  in  Reference  13.  The  simulation  studies  indicate  that  this 
level  of  acceleration  is  quite  satisfactory  and  therefore,  considerably 
smaller  winglets  could  be  employed  for  maneuvering  in  hover.  The  winglet 
required  for  level  trim  in  forv/ard  flight  with  a spherical  centerbody 
appears  to  result  in  larger  size  winglets  than  are  necessary  for  maneu- 
vering around  hover. 
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TABLE  IV 

HOVERING  STABILITY  AND  CONTROL  DERIVATIVES  FOR  TANDEM, 
SINGLE  ROTOR  HELICOPTERS  AND  FULL  SCALE  AEROCRANE 

TANDEM  HELICOPTER  (Ref.  k) 


design  gross  iveight 

25T 

lOOT 

l/Ft-Sec 

O.OOU9 

0.0027 

M l/Sec 

q 

-0.11 

-0.89 

Mg  l/Sec^-In 

0.27 

0.16 

t 

L l/Ft-Sec 

-0.0076 

-0.0054 

L l/Sec 

P 

-0.57 

-0.92 

Lg  l/Sec  ®-In 

0.54 

0.39 

NOTE: 

Tandem  has  no 
significant  Pitch- Roll 
coupling. 


SINGLE  ROTOR  (Ref.  10) 


NOTE: 

Speed  stability  and 
dihedral  effect  not 
given  in  Ref.  10. 
Coupling  derivatives 

(M  , L ) estimated. 

P q 


^ Estimated  assuming  = 4.5,  Y = 

^Estimated  using  scaling  laws  of  Ref.  I ~ w^*^;  I ~ w^’^,  R ~ w*^'^ 


DESIGN  GROSS  ^VEIGHT 

20T 

(2) 

lOOT^  ^ 

M 1/Sec 

q 

-0.31^^^ 

-0.36 

l/Sec^-In 

0.21^^^ 

0.11 

M l/Sec 

P 

0.21^^^ 

0.25 

L l/Sec 

-1.4 

-2.7 

L.  l/Sec^-In 
g 

0.38 

0.32 

L l/Sec 

q 

-0.96^^^ 

-1.86 
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TABLE  IV  Continued 
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Figure  l6:  Comparison  of  Damping  and  Control  Sensitivity  of  Single 
Rotor,  Tandem  and  aEROCRAIE. 


1.0 


Figure  17:  Comparison  of  aERCCRANE  Camping  and  Control  Sensitivity 
with  Handling  Qualities  Data  and  MIL-H-85OIA. 
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CONCLUSIONS 


The  aEROCRANE  configurations  examined  in  this  report  with  one  minor 
exception  (with  the  sling  load  attached  at  the  center  of  gravity) 
are  dynamically  stable  in  hovering  flight.  The  sling  load  configur- 
ation has  only  a small  influence  on  the  dynamic  stability. 

Compared  to  large  proposed  helicopter  designs,  the  AEROCRANE  config- 
uration studied  has  very  high  levels  of  angular  damping  and  low 
control  power.  The  high  levels  of  angular  damping  result  from  the 
rigid  rotor  configuration  and  the  low  control  power  results  from 
low  rotor  angular  velocity.  The  pitch/roll  coupling  of  the  AEROCRANE 
is  similar  to  a single  rotor  helicopter. 

The  AEROCRANE  configurations  investigated  meet  the  visual  flight 
handling  qualities  specifications  of  MIL  H-85OI-A  in  terms  of 
damping  and  control  sensitivity.  It  should  be  noted  however  that 
application  of  MIL  H-85OI-A  to  large  AEROCRANES  requires  extrapo- 
lation to  gross  weights  well  beyond  the  region  where  experimental 
flight  test  data  are  available. 

A simulation  study  of  the  precision  control  characteristics  of  the 
AEROCRANE  in  hovering  flight  indicates  that  with  cyclic  control, 
the  task  of  hovering  the  AEROCRANE  over  a spot  in  the  presence  of 
gusts  or  translating  from  point  to  point  is  difficult.  The 
results  obtained  were  similar  in  many  respects  to  results  obtained 
in  simulation  studies  of  conventional  large  helicopters. 
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Somewhat  larger  displacements  were  noted  in  the  study  reported  here 
for  these  tasks  in  contrast  to  the  conventional  large  helicopter. 
However,  in  both  cases  the  displacements  were  comparatively  large. 

It  should  be  noted  that  a simple  simulation  was  employed  here  and 
the  operator  controlled  the  vehicle  from  a presentation  of  trans- 
lational displacements  only. 

5. )  The  simulation  study  also  examined  the  precision  control  character- 

istics of  the  aEROCRAHE  with  winglets  which  provide  direct  transla- 
tional force  control  in  hovering  flight.  Both  of  the  tasks,  hovering 
over  a spot  in  the  presence  of  gusts  and  translating  from  point  to 
point  become  quite  simple  to  perform.  It  was  possible  to  hover 
over  a spot  in  the  presence  of  5 fps  RMS  gusts  within  a one  foot 
square,  using  the  winglets  with  the  basic  vehicle  (no  feedback). 

This  performance  approaches  that  obtained  in  the  Heavy  Lift  Heli- 
copter studies  with  a very  sophisticated  automatic  control  system. 

6. )  Comparatively  large  winglets  are  required  to  trim  the  aEROCRANE  with 

a spherical  centerbody  in  level  flight  at  an  advance  ratio  of  0.3* 

The  simulation  studies  indicated  that  considerably  smaller  winglets 
could  be  employed  if  the  winglets  are  employed  for  maneuvering  in 
hover  only. 

7. )  Sling  load  dynamics  must  be  carefully  considered  in  the  control  of 

any  hovering  vehicle  which  obtains  an  appreciable  portion  of  its 
lift  from  buoyant  force  and  carries  a sling  load  which  is  an 
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appreciable  fraction  of  the  gross  weight.  This  is  particularly  true 
when  hovering  control  is  obtained  by  applying  moments  to  the  vehicle 
(e.g.  cyclic  pitch).  The  presence  of  the  buoyant  force  reduces  the 
translational  acceleration  obtained  from  pitch  attitude.  There  is  a 
range  of  control  input  frequencies  in  the  vicinity  of  the  sling  load 
frequency  where  as  a result  of  the  sling  load  motions  the  translational 
acceleration  arising  from  pitch  attitude  reverses  sign.  That  is,  a 
nose  doTO  attitude  produces  a rearward  acceleration.  There  is  strong 
coupling  between  the  translational  acceleration  and  sling  load  motion. 
The  effect  of  the  sling  load  may  be  approximately  treated  as  a forced 
response  problem  rather  than  dynamically  coupled  to  the  vehicle. 
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APPENDIX  A 


WINGLET  CONTRIBUTIONS  TO  CONTROL  AND  STABILITY  OF  THE  AEROCRANE 

An  alternate  means  of  controlling  the  AEROCRANE  has  been  proposed 
which  involves  the  use  of  winglets.  Winglets  are  vertical  surfaces 
mounted  at  the  tips  of  each  rotor  blade  as  illustrated  in  Figure  A-1. 

The  winglets  are  of  interest  for  two  reasons.  First  it  is  possible 
to  fly  forward  producing  a net  translational  force  by  suitably  varying 
the  winglet  deflection  angle,  , such  that  forward  flight  can  be 
achieved  without  tilting  the  vehicle.  Second,  the  winglet  can  be 
employed  for  maneuvering  the  vehicle  as  well,  by  providing  transla- 
tional forces  directly  through  winglet  deflection  rather  than  through 
the  tilting  of  the  thrust  vector  with  cyclic  pitch.  This  method  of 
control  looks  particularly  attractive  for  precision  hovering  control 
as  discussed  elsewhere  in  this  report.  Reference  ik  has  examined  the  use 
of  winglets  for  control  of  the  AEROCRANE  and  indicates  the  desirable  control 
characteristics  which  can  be  obtained  with  winglets. 

The  contributions  of  the  winglets  to  the  stability  and  control 
characteristics  of  the  AEROCRANE  are  developed  in  this  Appendix.  From 
Figure  A-1  the  lift  and  the  drag  forces  on  a winglet  can  be  expressed 
as 


“ "2  p ^ an (6^  + e) 


(A-1) 


The  resultant  forces  applied  to  the  vehicle  by  a single  winglet  in  the 
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rotor  axis  system  is  from  Figiire  A-1, 


(A-2: 


sin  (e  + t ) - Dw  (e  + \)r ) 

= - L„  cos  (e  + ^)  + sin  (e  + iji ) 

The  angle  between  the  resultant  velocity  experienced  by  the  winglet, 
V,  , and  the  velocity  due  to  rotation,  OR,  is  given  by 


u.  cos  llf 

tan  e = = ^ 

1 + p,  cos  l)f 


or 


li  . , 1 + u sin  il; 

sin  e = cos  it  and  cos  s = 

Ur  Ur 


(A-3) 


These  relationships  can  be  used  to  give 


. \ cos  t 

sin  (\j(  + e ) = ^ 


Ur 


\ Sin  + u 
cos  (;(r  + e ) = ^ ^ 

U R 


and 


U R ~ 


1 + u sin 
cos  € 


(A-U) 


The  velocities  have  been  nondimensionalized  by  the  tip  speed,  QR, 
such  that. 


U 


X 

QR 


and  p n 


Vr 

Sr  * 


Since  the  primary  interest  is  to  develop  the  contributions  of  the  winglet 
to  stability  and  control,  the  winglet  drag  force  will  be  neglected  in  the 
succeeding  development,  as  it  will  not  be  a primary  contributor  to  the 
forces  and  moments  influencing  the  stability  and  control  characteristics. 
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It  is,  of  course,  important  in  estimating  the  power  consumption  of  the 
winglets  to  include  the  winglet  drag  force. 

Neglecting  the  winglet  drag  and  introducing  equations  (a-1)  and  (a-^) 
into  (a-2),  the  resultant  forces  applied  to  the  vehicle  can  be  expressed 
as 


— 2p  ■*'s)  COS  \{r 


(A-5) 


Ykw  = - i P Srt  a„  hR(6  w + e ) (sin  t + d ) 

It  is  further  assumed  that  e is  a small  angle.  This  is  a very  satisfactory 
assumption  particularly  for  studies  near  hovering  flight,  and  probably  up 
to  an  advance  ratio  p,  = 0.3,  where  the  maximum  value  of  s is  23°  as  given 
by  equation  (A-3)«  It  is  further  assumed  that  the  angular  control  motion 
of  the  winglet  is  sinusoidal  such  that  the  angle  6^  is  given  by 


6,^  = 6^  cos  + 6^  sin  i|r  (a-6) 

With  the  assumption  that  e is  small  and  using  equation  (a-6),  equation 
(A-5)  can  be  expressed  as 

I 6 + p,  6 p, 

Xr«  = - i p (nR)^  S*,  aJ  ( ° g - ) + -IP  sin  t 

6 u + P- 

+ -jp  cos  \|r  + ( 2 — ) ^ 

6 u 6 a ) 

+ -p-  sin  3t p-  cos  3'>  > (A-7) 
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= - i D (nR)^  6^,  a,, 


6.(1  + 


+ (m,^  f 6^p,)  cos  \|t 


+ 6^|j,  sin  \|r  + i-  (m>  "*■  6^(1  + ))  sin  2ijr 


- i 6^|j,  cos  3t  " 4 63M‘  sin  3^  | 

Equations  (a-7)  give  the  resultant  forces  produced  by  a single  winglet. 

It  can  be  seen  that  control  deflections  6^  and  5^  produce  constant  forces 
as  well  as  fluctuating  forces  up  to  the  third  harmonic. 

If  four  winglets  are  employed,  the  resultant  forces  applied  to  the 
vehicle  can  be  calculated  by  summing  the  contributions  of  the  four  winglets 
using  equations  (A-7)  for  each  winglet  with  the  azimuth  angle,  \|( , suitably 
adjusted  to  account  for  the  azimuth  location  of  each  winglet,  i.e., 

ta  = I"!  90° 

if3  = + 180° 

'|t4  = Vi  + 270° 

All  of  the  harmonic  components  cancel  leaving  for  four  winglets 
^Rw  " ” P ^w(6(,  M* ) 

= - p (nR)=  a„(6  ) (1  + ) 


(A-8) 


(A-9) 


It  should  be  noted  that  while  the  harmonics  cancel  in  an  overall  sense, 
the  winglets  individually  do  apply  oscillating  forces  to  the  individual 
blades  as  given  by  equations  (a-7) 
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Thus  equations  (a-9)  give  the  forces  applied  to  the  AEROCRANE  as  a 
function  of  control  deflection  and  advance  ratio.  It  should  be  noted 
that  the  winglets  will  produce  moments  as  well  as  forces  since  the  center 
of  gravity  of  the  vehicle  is  in  general  not  located  in  the  plane  of  the 
rotor. 

The  size  of  the  winglet  can  be  related  to  the  requirement  to  trim  the 
vehicle  in  forward  flight  in  level  attitude  equilibrium.  Thus,  neglecting 
the  rotor  in  plane  force,  the  longitudinal  winglet  force  can  be  equated  to 
the  drag  of  the  centerbody  to  give  an  estimate  of  winglet  size  required 
- p (nR)®  S„  a„  (6^  + M.)  - 2 p n Rg^  0^=0 

Solving  for  winglet  deflection, 

6^  = - - 2 s (A-10) 

where 


TT  Rg^  Cg 


Taking  the  winglet  to  be  of  rectangular  planform  with  a chord  equal  to 
the  blade  chord,  c 


a a b^ 


where  b„  is  the  winglet  span  nondimens ionali zed  by  rotor  radius. 


In  a similar  fashion  the  sine  component  of  winglet  deflection  must 
balance  the  i«Iagnus  force.  Equating  the  side  force  produced  by  the  winglet 
to  the  Magnus  force  gives 


and  therefore,  the  relationship  between  the  two  components  of  winglet 
deflection  is  given  by  equation  (a-10)  and  (a-11)  as 


(1  + 


(A-12) 


The  area  of  the  winglet  should  be  selected  such  that  winglet  stall  is  not 
encountered.  The  winglet  angle-of-attack  is  given  by 


= e + 


(6  cos  i|r  + 6 sin  ilr 
c s 


+ 


u cos  'if V 

1 + p,  sin  f ' 


(A-13) 


Substituting  the  required  winglet  control  deflections  required  for  trim, 
equation  (a-13)  can  be  expressed  as 


- M- 


sin  t cos  „ 

y— — ^ r + 2s 

1 + p,  sin 


cos 


U + 0) 


where 


0 = tan 




1 + p® 


(A-lU) 


Figure  A-2  shows  the  winglet  angle-of-attack  variation  with  azimuth 
required  for  level  attitude  trim  at  an  advance  ratio  of  0.3.  The  sample 
calculation  is  based  on  the  following  dimensionless  parameters . 


Co  = 0.80 

a = .138 

^ U M ~ 30 

X = -^6 

a^  = 4 

This  vehicle  geometry  is  based  on  the  configuration  described  in 
Reference  1.  The  dreg  and  '-lagnus  force  coefficients  were  determined 
experimentally  as  described  in  Reference  1.  The  cyclic  amplitudes  of 
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winglet  motion  required  are  also  shown. 


The  analysis  assumes  that  the  winglet  does  not  stall  and  so  the 
results  of  Figure  A-2  indicate  that  the  winglet  span  to  radius  ratio, 
must  be  larger  than  0.25  to  maintain  the  winglet  angle-of -attack 
below  15°  and  consequently  avoid  stall.  If  the  winglet  is  also  used  for 
control  at  forward  speed  then  some  additional  margin  must  be  available 
for  control  in  addition  to  trim.  A further  characteristic  of  interest 
is  the  acceleration  capability  of  the  winglet,  particularly  with  respect 
to  maneuvering  near  hover.  The  ratio  of  the  winglet  force  produced  by 
winglet  deflection  to  the  buoyant  force  can  be  expressed  as 


, aXp„  3o-  a bj^ 

-L  2- 


Fa  96^ 


« ""c  16  g x" 


(Q^R) 


(a-15) 


Using  the  dimensionless  parameters  given  above,  and  expressing  (a-15) 
in  terms  of  vehicle  weight  and  per  degree  of  winglet  deflection 


W 


dX,  w 
d 6 


Fg 


o 1700  w 


(A-15 


Equation  (a-15)  is  an  expression  for  the  acceleration  in  g's  produced  by 
one-degree  of  winglet  deflection.  Linear  acceleration  is  independent  of 
size,  based  on  Froude  scaling.  Using  the  parameters  of  the  dynamic  model 
described  in  Reference  1 for  the  rotor  rotational  speed  and  radius, 


n = 30  RFM 

R = 19.9  Ft. 

equation  (A-15)  becomes 

1 3X,„  ^ Fg 

rr  _ = 0.12  b„  — 


(A-15) 
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For  a buoyant  force  to  weight  ratio  0.75,  typical  of  proposed  AEROCRAIIE 
designs  with  payload  and  a = 0.25  the  minimum  size  winglet  for  trim  at 
an  advance  ratio  of  0.3, 

f7  ^ = 0.023  g per  degree  Ca-16) 

^ S 6 ° 
c 

The  significance  of  this  level  in  controlling  the  aEROCRANE  near  hover 
is  discussed  elsewhere  in  this  report. 

Examination  of  equations  (A-9)  indicates  that  the  addition  of  winglets 
will  also  influence  the  stability  and  gust  sensitivity  of  the  AEROCRAME 
from  the  fact  that  the  longitudinal  force  depends  linearly  on  the  advance 
ratio.  In  hovering  flight  the  linearized  dependence  of  the  side  force 
on  advance  ratio  is  zero.  Note  that  the  linearized  contributions  of  the 
centerbody  are  also  zero  so  the  impact  of  winglets  must  be  carefully  examined 
not  only  with  respect  to  control  but  with  respect  to  gust  sensitivity. 

Consequently,  the  following  terms  must  be  added  to  the  equations  of 
motion  of  the  AEROCRANE  to  account  for  the  presence  of  winglets  in  near 
hovering  flight. 

Equations  (a-9)  are  given  in  wind  axes.  Expressing  the  contributions 
of  the  winglet  in  body  axes 

AXb«  = ^ ^ u 


(A-17) 


= rr-^s  ^ ~ ^ 

s 

where  6 and  6 are  now  assumed  to  be  body  referenced.  The  derivatives 
c s 

are  equal  to 
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(a-18) 


3 6 3 6 

c s 


3X, 


3 V 


p (QR)  S a 


In  addition  to  forces,  the  winglets  also  produce  moments  owing  to  the  fact 
that  the  center  of  gravity  of  the  vehicle  is  located  a distance  r^  below  the 
rotor  plane.  The  pitching  moment  and  rolling  moment  produced  are 


ALe  K - r 
1 ° 


(a-19^ 


Thus  control  moments  and  contributions  to  the  speed  stability  and  dihedral 
effect  are  produced.  Using  equations  (a-17),  (A-18)  and  (A-19)  these  are 

3AM0H 

— = p(nR)^  S„  r^  a„ 


= p(QR)  r^  a.  ^ 


(A-20) 


3AL3  ^ 
3 6_ 


3AL9H 

3 V 


-p(r2R)^  Sm  r^  aw 


-p(nR)  Sw  aw 


The  influence  of  these  various  winglet  contributions  is  discussed  in  the 
main  body  of  this  report. 


llU 


OP 

<^T7^CK. 

Poe.TU.j(v\^ 

DE<^ 


\^I^^K,V^  CVOUC.  *^rrtPL\’TUDE  PoO.  UE'ieU 

AiTX  ITUDe"  T^irA  t 0 i O , <^=o) 


Figure  A-2.  Winglet  Angle-of-Attack  Required  for  Level  Attitude 
Trim  (9  = 0,  0_=  0)  as  a Function  of  Winglet  Span  to 
Rotor  Radius  (bj^)  at  an  Advance  Ratio  of  u = 0.3. 
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APPENDIX  B 


ROTOR  AERODYNAMIC  MODEL 

The  aerodynamic  forces  and  moments  produced  by  the  rotor  are 
modelled  based  on  the  following  assumptions  in  addition  to  the  usual 
assumptions  that  the  blade  inflow  angle  is  small  and  the  dependence 
of  blade  element  lift  on  angle-of-attack  is  linear 

1. )  The  rotor  blades  are  assumed  to  be  rigid  so  that  there  is 

no  flapping  motion.  This  assumption  is  a consequence  of  the  method  of 
blade  support  as  well  as  the  low  angular  velocity  of  the  rotor.  It  is 
possible  that  for  a very  large  AEROCRAHE  this  assumption  may  not  be 
valid.  However,  it  is  necessary  to  have  detailed  designs  for  the 
blades  and  blade  support  system  before  the  validity  of  this  assumption 
can  be  examined  further  as  a function  of  vehicle  size. 

2. )  The  rotor  induced  velocity  is  modelled  in  the  following 
way.  The  induced  velocity  is  assumed  to  consist  of  component  inde- 
pendent of  azimuth  and  first  harmonic  cos  Lie  and  sine  components 
such  that  the  total  inflow  velocity  can  be  written  as 

\ = Xq  + (XiX  + X m)  cos  i|f  ^ l sin  ijr 

The  azimuth  independent  component,  X^,  is  assumed  to  be  given  by 
momentum  theory  based  on  an  actuator  disc  area  which  does  not  include 
the  spherical  centerbody.  This  component  is  assumed  to  be  independent 
of  radius.  Although  this  assumption  is  inconsistant  with  the  fact 
that  the  blades  are  untwisted  (Ref.  5)  it  should  be  satisfactory  for 
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prediction  of  stability  and  control  characteristics . A more  refined 
theory  is  desirable  when  predicting  performance,  i.e.,  power  required. 
Thus, 


\ = 
0 


The  harmonic  components  of  the  induced  velocity  arise  from  two  sources: 
the  blow  back  of  the  rotor  wake  which  results  from  vehicle  translation 
and  results  in  the  component  X ^ x cos  \jf ; and  the  self  induced  components 
X M and  X ^ which  arise  from  the  fact  that  the  non-flapping  rotor  produces 
aerodynamic  hub  moments. 

Xi,  the  blow  back  effects  is  assumed  to  be  given  by  doubling  the 
value  obtained  from  Coleman’s  theory.  That 


Xi  = - 2 X^  tan  | 


where 


X = tan 


The  doubling  of  the  theoretical  value  is  based  on  the  correlation  of 
experimental  results  with  theory  presented  in  Reference  1 as  well  as 
other  indications  from  the  literature. 

X M aiad  X i.  are  calculated  by  assuming  that  they  are  proportional 
to  the  aerodynamic  hub  moments  produced  by  the  rotor.  That  is, 

2C„ 

2Cu 
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J 


The  proportionality  factor  between  the  aerodynamic  hub  moments  and 


the  harmonic  inflow  components  was  assumed  to  be  given  by 


. _ X /aa  N 

' ■ 2 ‘ra  ’ 


This  result  can  be  derived  from  momentum  theory  as  shown  for  example  in 
Reference  8.  The  constant  in  this  expression  depends  upon  the  assumptions 
employed  in  the  momentum  analysis.  The  value  of  i for  the  constant  in 
this  expression  was  initially  based  on  the  experimental  results  of  Reference 
7.  It  is  shown  in  Reference  1 that  this  model  gave  good  agreement  between 
the  measured  and  predicted  dynamic  response  characteristics  of  the  AEROCRAKE. 
Reference  9 also  presents  experimentally  determined  values  of  this  quantity 
based  on  frequency  response  tests  of  a model  rotor.  The  results  of  Reference 
9,  given  in  terms  of  the  parameter  L,  where  j = Ly , agree  well  with  the 
above  expression  for  j. 

The  only  other  feature  of  the  rotor  force  and  moment  prediction  to 
be  noted  is  the  fact  that  the  rotor  blades  end  at  the  radius  of  the  spherical 
centerbody  (Rg)  and  consequently  various  factors  appear  in  the  rotor  aero- 
dynamic equations  to  account  for  this  fact.  Define 


The  equations  for  the  rotor  forces  and  moments  are 


^s  ■ '^s 


p ^ u M-o  ^1 
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Xl  fj 


IJ,  X 
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■ ^ 

Xi  f4  H 
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2C, 
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“5“  [^4  + 2 M-< 


3]  + iJ-, 


The  dominant  rotor  contributions  to  the  vehicle  stability  and  control 

2C  „ 2C , 

arise  from  the  hub  moment  contributions  given  by  and  . It  can 

° a CT  a a 

be  noted  from  the  above  expressions  that  a number  of  the  rotor  hub 
moment  derivatives  are  similar  in  hovering  flight.  If  the  blow  back  effect 
is  approximated  at  low  speeds  by  assuming  the  wake  angle  small 
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and 


X 1 = - U, 


As  a consequence  of  this  assumption  it  can  be  seen  from  the  above 
expression  for  the  pitching  moment  that 


-2. 

a a 3Bis 


_2_ 

^ _ _2_ 

_2_ 

_2_ 

3Cl 

aa 

3 u 

aa 

Bq 

aa 

BAis 

aa 

Bp 

_2_ 

aa 

3Ci. 

3 u, 

can  be 

expressed 

in  terms 

of 

the  ■ 

as  follows.  In  hovering  flight 


2Ct 

acr 


^■3 


and 


_2_ 

aa 


BCl 
3 u 


_ 3 


^3  + fj 


1 + 


Substituting 


_2_ 

aa 


3Cl 
3 u. 


2C  T X 
aa  ^ 

l + f f3 


Taking  j = ^ selecting  values  for  solidity  and  ratio  of 

I^ol 

centerbody  radius  to  rotor  radius  typical  of  AEROCRANE  designs,  the 
hub  moment  derivatives  can  be  calculated  as  a function  of  thrust 
coefficient.  These  terms  are  given  in  Figure  B-1. 
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APPENDIX  C 


AEROCRAm  EQUATIONS  OF  MOTION  WITH  SLING  LOAD 

Reference  1 presents  a detailed  derivation  of  the  AERCCRANE  equations 
of  motion  with  a sling  load.  The  sling  load  is  taken  to  be  supported  at 
a distance  Z^  below  the  center  of  gravity  of  the  vehicle.  The  sling  load 
is  assumed  to  be  a point  mass  with  negligible  aerodynamic  forces  acting 
on  it  supported  on  a cable  of  length  Z,_. 

The  six  equations  of  motion  describing  the  four  degree-of-freedcm 
dynamic  motion  of  the  AEROCRANE  in  hover  with  a two  degree-of-freedom 
sling  load  are 

Fb  0 - Za  m,.  g (35,.  = Ls/;, 

Fb  9 - Za  m,.  g = Mj/a 

(0  (Za  + Zj  + 'i  Zj  = 0 

a (Za  + Zj  + 0 Z J = 0 

9l  + (Wq  - Fg)  9 = Xg/A- 

0L  - (Wq  - Fg)  0 = Yg/A 

The  terms  represented  explicitly  above  include  the  inertial  and 
gyroscopic  terms,  buoyancy  and  gravity  terms  and  aerodynamic  acceleration 
forces  and  moments  ("apparent  mass"  terms),  m^  also  includes  the  apparent 


I^  9 + I n 9 + r mA  V + 
z 0 ' s 

-I  n0  + I^9-r  hiaU  + 
z o ^ s 

- m^  Ug  - m^  g (9  + 9l  ) - 

- m^  Vg  + m,.  g (0  + 0L  ) + m,. 

m'Ug  - ffiA  *9*  - m^  g 

m^v  + r mA0+m,.g 
s o 
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mass  of  the  centerbody  and  includes  its  effect  cn  the  pitching  and 
rolling  inertia.  The  right  hand  side  terms  include  all  of  the  aerodynamic 
forces  and  moments.  The  axes  and  geometry  are  shewn  in  Figure  C-1. 

The  aerodynamic  forces  and  moments  arise  from  the  rotor  center- 
body  and  winglets. 

The  moments  can  be  expressed  as 

Lj/*  = Lh  + r^  (Y,  + + F„  ) 

M,/*  = Mh  + r^  (H,  - + D) 

The  forces  are 

Xs/*  = - H,  + - D 

Ys/a  = Y,  + Y^  + F„ 

where 

Lh;  Mh  are  rotor  hub  moments 
Yf,,  H,  are  rotor  inplane  forces 

X^,  Y^  are  winglet  forces 

Fm , D are  centerbody  forces,  that  is,  Magnus  force  and  drag. 
Appendix  B presents  expressions  for  the  rotor  forces  which  can  be  used  to 
determine  the  rotor  derivatives  and  Appendix  A presents  the  expressions 
for  the  winglet  contributions.  The  drag  and  Magnus  forces  are  given  by 

D = i p nRf  C,  (u^  + v^) 

~ 'k  Q ■n'Rf  C|_ni  (u'  + v^) 
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where  Cq  and  C,  were  determined  in  Reference  1 by  correlation  with 

M 

experiment.  Note  that  if  these  terms  are  linearized  then  their  contri- 
bution to  the  linear  equations  of  motion  are  zero.  In  the  analog  computer 
study  discussed  elsewhere,  these  non-linear  terms  are  retained  in  the 
computation. 

The  equations  of  motion  are  now  non  dimensionalized  using  the 
approach  described  elsewhere  in  this  report. 

The  first  two  equations,  the  vehicle  moment  equations,  are  divided 
by  Q,  and  is  included  in  the  acceleration  terms  to  nondimensionalize 
the  time.  The  moment  equations  become 
I _ 

0"  + ^ 9'  + p r^  (2,2  -/  + 9 (2^3  0 ■_  y,  (11^2  = aCl 

I 

- P 0'  + 9"  - {i  0 i-o  u'  + P (2  9 - li  Z,  (S  ,2  9 ^ = AC  „ 


The  vehicle  force  equations  are  divided  through  by  m^  R to  give 


2 ^o^ 


B"  - 


^ 5-  9 L + 9 = B C. 

F,  (1  + |)  F«  (1  + |) 


V + 
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S - 

2 ^0  ^ 
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0"  + 


Fr  (1  + |) 


0L ~ ■ 0 = B C 

Fr  (1+1)  ^ 


where 


A k 2 

T,  _ y 2 

rrs 


The  sling  load  equations  are  divided  through  by  the  quantity  m,_g.  The 
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signs  are  reversed  in  the  first  sling  load  equation 


The  first  grouping  of  terms  are  the  rotor  derivatives  given  in  Appendix  B 
The  winglet  terms  are  given  in  Appendix  A,  and  they  are  expressed  nondimen 
sionally  near  hover  as 


S a 
_ w w 

■ tt  R = 


(6^  + M.) 


S a 
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As  a result  of  symmetry  there  are  a number  of  equalities  among  the  rotor 
derivatives  as  indicated  below  in  the  list  of  derivatives  which  are  found 


to  be  present  from  the  results  presented  in  Appendix  B. 
Cl.  = C y 

C^.  =(!,_=-  C„ 


c 


c 


H 


a 


IS 


The  winglet  contributions  are  given  above  where  in  hovering  it  may  be 
noted  that 


and 
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The  complete  equations  of  motion  for  the  hovering  aEROCRANE  may 


be  written  in  matrix  notation  as 

[M]  Cq"]  + [C]  [q'}  + K {q}  = F [6} 
where  the  motion  variables  are 


and  the  control  terms  are 


The  matrix  elements  for  the  dimensionless  equations  are  given  on  the 
following  pages. 
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